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Abstract 

It is well known that an elastic sheet loaded in tension will wrinkle and that the length scale of the 
, wrinkles tends to zero with vanishing thickness of the sheet [Cerda and Mahadevan, Phys. Rev. Lett. 

CLf 90, 074302 (2003)]. We give the first mathematically rigorous analysis of such a problem. Since our 

methods require an explicit understanding of the underlying (convex) relaxed problem, we focus on the 
wrinkling of an annular sheet loaded in the radial direction [Davidovitch et al., PNAS 108 (2011), no. 
45]. Our main achievement is identification of the scaling law of the minimum energy as the thickness 
of the sheet tends to zero. This requires proving an upper bound and a lower bound that scale the 
same way. We prove both bounds first in a simplified Kirchhoff-Love setting and then in the nonlinear 
three-dimensional setting. To obtain the optimal upper bound, we need to adjust a naive construction 
(one family of wrinkles superimposed on a planar deformation) by introducing a cascade of wrinkles. 
The lower bound is more subtle, since it must be ansatz-free. 

1 Introduction 

In the last few years the wrinkling and folding of thin elastic sheets has attracted a lot of attention in both 
the mathematics and physics communities (see, e.g., the recent book by Audoly and Pomeau Wrinkled 
configurations can be viewed as (local) minimizers of a suitable elastic energy, consisting of a non-convex 
"membrane energy" plus a higher-order singular perturbation representing "bending energy" . Though the 
physically relevant wrinkled configurations are local minimizers, we can begin to understand their character 
by focusing on (i) the minimum value of the elastic energy, and (ii) the properties of low-energy deformations. 
In this paper we identify the scaling law of the minimum energy for an annular sheet stretched in the radial 
direction. This requires proving an upper bound and a lower bound that scale the same way. A naive 
approach to the upper bound, based on a single length scale of wrinkling, fails to achieve the optimal 
scaling [TU]; the successful approach uses a cascade of wrinkles. The lower bound is more subtle, since it 
must be ansatz-free. We prove it first in a reduced Kirchhoff-Love setting and later in a general ansatz-free 
three-dimensional setting. 

As mentioned above, the behavior of thin elastic sheets has attracted considerable attention from the 
physics community (see, e.g., work on sheets of graphene [H]). Mahadevan and Cerda considered the 
stretching of a rectangular elastic sheet with clamped boundaries [S] (see also [3S] for experiments and [T3] 
for numerical computation), by minimizing the elastic energy of the membrane within a particular ansatz. 
The problem we consider here is similar, but our viewpoint and achievement are different: we prove an upper 
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bound and a matching ansatz-free lower bound on the elastic energy. Our analysis does not assume a specific 
form for the solution. 

Our treatment requires knowledge of the underlying convex relaxed problem. In the radial setting 
(see [lOj ) the relaxed problem reduces to a simple one dimensional variational problem which can be ana- 
lyzed quite completely. For this reason we focus on an annular sheet stretched in the radial direction as a 
convenient model problem, rather than addressing the case considered in 

One might ask why we are so interested in the scaling law of the minimal energy. As mentioned above, 
stable configurations are local minimizers of the elastic energy. It seems difficult to find such configurations 
analytically (the associated fourth-order PDE is highly nonlinear). But we expect the physically-relevant 
configurations to have relatively low energy. Therefore we can obtain some information about them by 
identifying the energy scaling law, then investigating the properties of configurations that achieve this law. 
While the present paper focuses mainly on the energy scaling law, certain consequences are immediately 
evident. In particular, since our scaling law is linear in h and the bending energy is ft,^J|V^U3p, it is 
immediately evident that the low-energy configurations become increasingly complex as ft, — >■ 0. 

1.1 Context 

Motivated by experiments, several physics papers have studied the wrinkling of thin elastic films from a 
theoretical point of view. As already mentioned, Davidovitch et al. |10| considered an annular film stretched 
in the radial direction (see also (TS] for related results). Dead loads applied both on the inside and outside 
boundary cause the film to wrinkle in some region. Indeed, if the loads inside are large enough compared to 
the loads on the outer boundary, the deformation in the radial direction forces the concentric circles of the 
material near the inner boundary to decrease their length by more than is required by the Poisson ratio of 
the material. Therefore, the membrane needs to waste this excess in the circumference either by compression 
or by buckling out of plane, contributing to the energy with some amount which depends on h. In |10| they 
found an optimal solution (using energy minimization methods) within a particular ansatz and using a linear 
stress-strain law, obtaining conclusions about the extent of wrinkled region and the period and amplitude 
of wrinkles. In the present paper we consider the same problem using a nonlinear 3D model. We will prove 
an upper bound and a matching lower bound without assuming any ansatz. 

Our problem seems related to the experiment reported in (TB] . It consists of a circular thin elastic film 
placed on a liquid substrate with a droplet on top of the film. In this case, the capillary forces at the boundary 
stretch the film in the radial direction and the capillary forces from the drop force the film to wrinkle. This 
experiment was studied theoretically in j^7\ (see also [9]), though a lot of questions still remain open. We 
believe that our methods may also be useful in the study of this problem. 

The idea of proving an upper bound and a matching lower bound for the minimum of the energy has a 
long history; see e.g., the work of Kohn and Miiller on a model for martensitic phase transformation [50] 
(see also [8| for subsequent progress). As in our setting, the energy in |20| was composed of a nonconvex 
function of Du singularly perturbed by a higher order term. In the setting of |20j the minimizer develops a 
fine branching structure. Similar phenomena are seen in uniaxial ferromagnets and type I superconductors 
(see, e.g., [6], [7], and [21]). In all these settings there is a "relaxed problem", whose minimizers are the weak 
limits of optimal configurations as ft — > 0. The minimal energy for ft > is that of the relaxed problem plus 
a small correction that scales with ft. One difference here is the special character of the singular perturbation 
- bending energy rather than surface energy - which leads to creation of smooth wrinkles rather than walls. 

The main focus here is the scaling of wrinkles associated with ft > 0. This is different from mere 
identification of the extent of the wrinkled region, which can be done by studying the relaxed problem or 
using the tension field theory. There is a lot of literature on this application of tension field theory (see, e.g., 
the 1961 NASA report [26], [12], or a recent work on balloons ^). 

1.2 The main idea 

Before starting with rigorous arguments let us outline our result and the main ideas of its proof. We will work 
variationally, considering the sum of the elastic energy of the thin sheet and boundary terms representing 



2 



the work done by the loads. We first consider a simpUfied two-dimensional setting where the elastic energy is 
further split into a membrane and a bending term defined in terms of the midplane deformation u : — >■ K'^. 
The membrane term is written as the integral of a reduced 2D stored energy density obtained in a systematic 
way from the original 2>D stored energy density. As the bending term we choose the L? norm of the second 
derivatives of the out-of-plane displacement as often seen in the linear Foppl-von Karman energy used for 
small slopes and deformations. This 2D model is a curious hybrid since we use a nonlinear stretching term 
together with a linear bending term. One can ask why we don't also use a linear stretching part? In fact, 
that scenario would be very limited since it would lead to a very restrictive linear model for the relaxed 
problem. 

Since our focus is the limiting behavior as the thickness h of the sheet tends to zero, we divide the energy 
by h to get the energy per unit thickness Eh- The first step toward identification of the scaling law is to 
separate the contributions to Eh from wrinkling and from the bulk deformation. This is done by considering 
a relaxed problem, where instead of the original stored energy density we use its quasiconvexification and 
formally set ft, = (see [23] for more detail on this topic). 

Under mild assumptions the relaxed energy is convex, smaller than the original energy, and independent 
of the thickness. Moreover, we show that it possesses a unique solution uq (up to a translation) which is 
radially symmetric and planar. We denote the relaxed energy of uq by £q] as we will see the energy Eh{uo) 
is strictly larger than Sq- This is because uq involves compression in the hoop direction in a region close to 
the inner boundary (we will call it the "relaxed" region). Since the thickness is small the sheet prefers to 
wrinkle rather than to compress. 

The idea of the construction for the upper bound on the minimum energy is to superimpose wrinkles upon 
uq. After optimizing the amplitude and wavelength of the wrinkles in a naive ansatz we obtain a solution 
with energy Eq + Ch |log h\. To remove the logarithmic factor (i.e. to get the same scaling as the lower bound) 
we need to work harder. We observe that the out-of-plane part of the deformation decreases suddenly at the 
boundary of the relaxed region; this is the source of the | log h\ factor. At the same time, since the amplitude 
of the out-of-plane deformation is vanishing, the bending is vanishing as well. Therefore we can introduce 
branching of the wrinkles (changing their period) near the boundary of the relaxed region; this increases the 
bending but at the same time decreases the amplitude of the out-of-plane deformation, making the decrease 
of the out-of-plane displacement less steep. By arguing this way we obtain a construction whose energy is 
bounded by Sq + Ch. 

The lower bound miwEh > EQ + ch {c> Q independent of h) is proved using an argument by contradiction. 
In the simplified two-dimensional setting, we first use the relaxed problem to prove an estimate on the out- 
of-plane displacement M3. Then using interpolation we show smallness of Du^, allowing us to project the 
solution into the plane without changing its energy too much. Finally, we compare this projection with uq 
and obtain a contradiction from an argument about the area of the deformed annulus. The generalization 
to the nonlinear three-dimensional setting uses the main arguments of the 2D setting coupled with rigidity 
estimates first derived by Friesecke, James, and Miiller [T4] . 

The paper is organized as follows. In Section [5] we describe the three dimensional energy together with 
a reduced two dimensional model. The definition of the relaxed problem and a theorem about its unique 
minimizer is in Section [31 Section [31 contains both the upper bound and the lower bound in the reduced 2D 
setting. In Section [5l we generalize the upper and lower bounds proved in Section [31 to the three-dimensional 
setting. The last section contains a brief discussion of our achievements together with some open questions. 

We will use notation x — {xi,X2,xz) = [x' ^x^) for points in , A : B = tr(A^B) for the Frobenius inner 
product on matrices, and 9i, 82 for partial derivatives with respect to the first and second variable. 



We are interested in deformations of isotropic elastic thin films of annular shape. We consider a nonlinear 
three-dimensional elastic energy (per unit thickness) in a cylindrical domain with small thickness h 



2 The model 
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The stored energy density W^oiM) is assumed to be isotropic, so it can be written as a symmetric fmiction 
of the eigenvalues of VM^M : 

W3I5(M) = /3Z5(Ai,A2,A3). 

Here and below, we assume that M — Du has strictly positive determinant; this is natural, since M is the 
gradient of an elastic deformation. 

As already mentioned, we are interested in deformations of annular thin films. We consider a thin 
cylindrical domain Vt x (— ft,/2,/i/2) with a cross-section 

17 = {x e ]R2 . ^^^^ < ^^^^1 

for some radii < i?i„ < Rout ■ The dead loads are applied on the inner boundary in the radial direction (with 
magnitude T!i„, pointing inwards) and on the outer boundary (with magnitude To„t, pointing outwards), so 
the film will mostly stretch in the radial direction. These loads contribute to the total energy as 

where x — (a::i,X2,0) and dS denotes surface measure. We will show in Theorem [3] that under suitable 
assumptions on the elastic energy density Wj,d, radii Rin,Rout, and forces Tin, Tout, we have 

minEl^iu) + Bl^{u) = So + 0{h), (2.1) 

u 

where Eq is a constant (depending on Q,, Tin, Tout and W'^u)- Since is the limiting energy as /i — > 0, we 
view it as representing the "bulk energy" of the deformation. The order-/i correction is the contribution 
from the wrinkling of the membrane. 

2.1 The reduced model 

Since we consider domains which are thin in the zs-direction, we can gain insight by first considering 
a reduced two dimensional Kirchhoff-Love model. In this setting we are interested only in the deformation of 
one cross section (e.g. the mid-plane 0:3 = 0), knowing that we can extend the deformation to the thin three- 
dimensional body by assuming that straight lines normal to the plane remain straight and normal to the 
plane after deformation. Assuming this ansatz, the energy per unit thickness is the sum of the "membrane" 
and "bending" energies 

/ W{Du)dx' + h'^ j Q{Dv)Ax', 
Jn Jn 

where v is the normal to the mid-surface of the deformation, Q is a certain quadratic function (derived from 
Wsd) and the form of W{Du) will be discussed in Section [521 The second term (the bending energy) can 
be expressed using the first and second derivative of u. In the 2D analysis we will replace the bending term 
Q(Dv) by a simpler term |Z?^U3p. Though the new term doesn't represent a physically correct bending 
energy, it is mathematically more convenient and still captures the main phenomenon. After including the 
boundary terms the two-dimensional energy has the form: 

Eh{u):^ I W{Du) + h'^\D^U3fdx'+B{u), (2.2) 
Jn 

where the boundary terms are 

B{u) Trn f u{x') -^dS- Tout I u{x') ■ -f- dS. (2.3) 

J\x'\=B.i„ ^in J\x'\ = Raut tiout 

As in the general 31? setting our main result is a scaling law for the minimum of the energy. We will 
show that if Eh is defined by (|2.2p then 

min£;,,(u) =£o + 0{h) 

u 

for any sufficiently small h > 0. The constant £0 is the minimum of the relaxed energy (the same constant 
as in (j2.ip ). The order-ft, correction is the contribution from the wrinkling of the membrane. 
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2.2 The energy density 

In this section we will describe the assumptions we impose on the elastic energy density W^j). 

Since the energy density W^niM) is isotropic, it is convenient to represent it as a function of the principal 
strains (i.e., the eigenvalues Ai > A2 > A3 > of (M^Af)^/^): 

W^3d(M) = /3d(Ai,A2,A3). 

We assume that 

/313 e C2([0,00)3), /3^(1,1,1) = 0, /3D>0, 0<D^hD<C, 

/3i3(Ai, A2, A3) > Co(A? + A2 + Al)f/2 - Ci, 

where p £ (1, 2]. 

The main motivation to consider a p-th power lower bound for the energy density for large strains rather 
than the quadratic growth is to include a broader range of materials. For example, Agostiniani et al. [T] 
showed that the growth condition in (j2.4p with p = 3/2 is satisfied by both the neo-Hookean compressible 
model and Mooney-Rivlin compressible model, whereas these models do not satisfy quadratic growth for 
large matrices (i.e. (12. 4p with p = 2). 

To define an elastic energy in the reduced two-dimensional setting (i.e. for a map m : — > R'^), we 
need to define a stored energy density as a function of Du, i.e. for 3x2 matrices. One way to do this is to 
optimize the missing third component. We set 

W(M) := min WsoiMin, (2.5) 

where M S R'^^^ and M\^ denotes a 3 x 3 matrix with first two columns identical with M and ^ as the third 
column. It turns out that if we write 

J := det(F), C := F^F, h ■= ^'"^ tr(C), h {{^AC)? - tr(C2)) , 

and if 

|f(/i,/2,J)>0, |^(/i,/2,J)>0, |^(/i,/2,J) + |^(/i,/2,J) >0, (2.6) 

Oil oh oh Oh 

then the ^ that achieves the minimum in (|2.5p has to satisfy ^ _L M . For completeness, we give a proof of 
this fact in the Appendix (see Lemma I A. 

Isotropy of the energy density implies that W{M) is a symmetric function of eigenvalues of ^/ M'^M: 

W{M) = f{M,\2). (2.7) 

It is easy to see that the function / is related to /3D. Let 

w'3d(Ai,A2) := argmin /3£)(Ai, A2,i). 
t>o 

Then we immediately obtain 

VF(M) = /(Ai,A2) = /3d(Ai,A2,w;3I3(Ai,A2)). 
Moreover, the function / inherits properties of /3D. Indeed, (12. 4p implies that 

/eC2([0,oo)2), /(1,1) = 0, />0, 0<D'f<C, 
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In contrast with three dimensions, the case p = 2 in (j2.8p is not very restrictive in two dimensions. For 
example, two-dimensional energy densities obtained from incompressible three-dimensional models have often 
quadratic growth at infinity, and so they satisfy quadratic lower bound for large strains. We will prove our 
main results assuming p = 2 in two dimensions and 1 < p < 2 in three dimensions. 

For a given Ai > 1 we define w{Xi) as the point of minimum for the function /(Ai, •), i.e. 

f{Xuw{X,)) = mm f{X,,t); (2.9) 

we call w{Xi) the natural width of the strip with first principal strain Ai. We assume that for Ai > 1 

w{Xi) is a differentiable and non-increasing function. (2-10) 

We also assume that for Ai > 1 and A2 > w{Xi) the following conditions hold: 

(5i/(Ai, A2) - d2f{Xi, A2)) (Ai - A2) > 0, (2.11) 
dnf{Xi,w{Xi)) + di2f{Xi,w{Xi))w'iXi)>0, (2.12) 
ai2/(Ai,A2) >0. (2.13) 

The meaning of these relations will become apparent in a moment. Briefly stated we use them to show 
convexity of the relaxed energy (see Section |3]). The strict inequality in (|2.12p is not a typo — it is associated 
with strict convexity of the (relaxed) energy density in 2D in the tensile direction. 
Finally, we assume for Ai > 1 that 

detD^f{Xi,w{Xi))-(Xi-w(Xi)) > 9i2/(Ai,ii;(Ai)) •9i/(Ai,u;(Ai)). (2.14) 

Unlike (|2. 11112. 13| . this inequality does not seem to have a simple interpretation; however it is satisfied by 
typical choices of / (e.g. the one associated with an incompressible neoHookean 3D model). Condition (j2.14|) 
will be used in our analysis of the relaxed problem (Lemma 13.71 and equation (|3.17p ). 



3 The relaxed problem 

In this section we study relaxed problem and the properties of its minimizer. 

We define the relaxed energy density Wr{M) as the quasiconvexification of W{M) (see Pipkin [M] for 
more details) and the relaxed functional as 

Eq{u):= I Wr{Du)dx' + B{u), (3.1) 

where B{u) was defined in (|2.3p . 

It will be crucial to our analysis that Wr{M) is a convex function of the 3x2 matrix M. Pipkin proved 
in [23] that this is true whenever the unrelaxed density W{M) is a convex function of M. This is true 
for a broad range of material models. We shall assume throughout this paper that 

Wr{M) is a convex function of M. 

As with W^D and IF, it is convenient to represent the relaxed density Wr as a function of the principal 
strains: 

W^r(M) =/,(Ai,A2). 

We would like to write down explicitly. To do that we follow the idea of Pipkin j24j . Using the natural 
width w{X) defined by (|2.9p . we define 

/(Ai,A2) Ai > w(A2) and A2 > w(Ai), 

/(Ai,w(Ai)) A2<w(Ai) and Ai > 1, ^-j 
/(A2,w(A2)) Ai < w(A2) and A2 > 1, 
Ai < 1 and A2 < 1, 



/m(Al, A2 
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and WmiM) := /„i(Ai, A2). Pipkin showed in [57 that 

WriM) < Wm{M). 

We win show in a moment that under our hypotheses W„i{M) is convex, in particular Wm{M) < Wr{M), 
from which it follows immediately that Wm(M) = Wr{M). 

We want to show that Wm is convex. Pipkin [24] showed that this is equivalent to showing that the 
function /„ is convex and monotone in both variables: 

D^frn>0, <9<,/™>0, a = 1,2 (3.3) 

and satisfies the ordered force condition 

(9i/m(Ai,A2)-a2/™(Ai,A2))(Ai-A2) >0. (3.4) 

(When WsD is the energy density of the incompressible neo-Hookean material, f„i takes a particularly simple 
form, studied in [24g). 

It is natural to express (13.31) and (j3.4[) using some conditions on /. First, in the region where Ai > 
w{X2), A2 > w{Xi) condition (|3.3I) follows from (|2.8p . whereas the latter condition p.4[) is equivalent to (|2.11l) . 
In the second case of p.2p . we see from (12. 8p that 

9i/™(Ai,A2) - dif{\i,w{Xi))+d2f{Xi,w{Xi))w'{Xi)) = dif{Xi,w{Xi)) > 0, 

V ' 



and obviously 92/m(Ai, A2) = 0; therefore (13. 4p is satisfied in this case. By (|2.12l) we see that 

aii/™(Ai,A2) -aii/(Ai,u;(Ai)) + ai2/(Ai,«;(Ai))w'(A) >0, (3.5) 

and by definition of /,„ also 9i2/r(Ai, A2) = 922/r (Ai, A2) = 0. Since fm is symmetric, the third case follows 
immediately. 

Finally, in the last case Ai < 1, A2 < 1, both p.3p and (|3.4|) are trivially satisfied. We have shown that 
if / satisfies (|2.8|) . (|2.1ip . and ()2.12|) . then satisfies ()3.2|) . in particular Wm is convex. Therefore fr = fm 
and fr has the form p.2|) . 



3.1 The one-dimensional variational problem 

We want to find a minimizer of the relaxed energy p.ip . Assuming it is radially symmetric, we formulate a 
one-dimensional variational problem, which admits a unique minimizer v. Afterward we show some properties 
of V. 

To look for a radially symmetric minimizer, we consider 

uo{r,e)^{v{r),e) (3.6) 
in polar coordinates. Then p.l|) becomes the one-dimensional variational problem 



The function /r(Ai, A2) is defined for Ai > and A2 > 0. Since we do not assume a priori that > or 
w > a.e., we need to extend the domain of /,.. It is convenient to do it in the following way: 

^ ^^|/.(Ar,^.(A,)) ifA,>0,A2<0, 

^ ' \/.(|Ai|,A2) if Ai <0. ^ ' 

Under our assumptions both Wr and fr are convex functions, and so we expect this variational problem 
to be solvable using direct methods of Calculus of Variations provided TmRin < ToutRout- 
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Remark 3.1. If TinRm < ToutRout, we claim that any minimizer v of p.7p has tensile hoop stress some- 
where in [Rim Rout)- Indeed, if we denote by ar and ag radial and hoop stress, respectively, the optimality 
condition reads (rcr^)' = ae, ar{Rin) = Tin,Or{Rout) = Tout- Integrating the equation gives 



t-Raut 

I (Tg{r) dr = ToutRout - T^nRin > 0, 



and so (Je{r) > for some r G {Rm, Rout)- 

Remark 3.2. In the case TinRm = ToutRout, the Euler- Lagrange equation implies that the hoop stress is 
identically zero, and a minimizer of p.7p is unique only up to an additive constant (i.e. only v' is uniquely 
determined). 

If TinRin > ToutRout, H is casy to see that the energy in (|3.7|) is not bounded from below and so the 
minimization problem has no solution. 



Remark 3.3. It is important to understand what are the consequences of extension p.Sp . We observe that 
fr being even in Ai means that 9i/r(Ai, A2) < for Ai < 0. Therefore, 

9l/r(Al,A2) >0 

implies 

Ai > 0. 

Also, for any A2 < 0, the hoop stress 

ere = 92/r(Ai, A2) = 

no matter how large IA2I is. 

From now on we will always assume that 

TinRin ^ ToutRout- ('^■^) 

We claim that under this condition there exists a unique solution v to the variational problem p.7p . 

Theorem 1. Let f satisfy (j2. 8112. 14)) and let < Rin < Rout o,nd Tout be fixed. Then there exists a range 
of inner-boundary loads Tin, 0, subset of {Tout, Tout Rout /Rin), such that minimizer v of p.7p exists, it is 
unique, and the following holds: 



There exists L G {Rin, Rout) such that 



^ < w{v'{r)) r G {R^n,L), ^ > w{v'{r)) r G {L, Rout), 
r r 

i.e. there is tensile hoop stress in {L,Rout) whereas in {Rin,L) there is compression in the hoop 
direction. We will call the region with a tensile hoop stress a non-relaxed region and its complement a 
relaxed region. 

• The deformation v avoids interpenetration, i.e. 

v{r) > 0, v'{r) >Oforre {R^n,Rout)■ 

• Consider the function h{r) — w{v'{r)) — representing the amount of arclength we need to waste 
in the relaxed region (_/?,„, L). Then 

h'{L) < (3.10) 

(and obviously h{L) — 0.) 



8 



Remark. Condition p.lOp means that the excess arclength associated with wrinkling at radius r grows 
linearly as a function of the distance from L; we will see later that this introduces some difficulties in the 
upper hound. 

Remark. It is easy to show that the relaxed energy associated with an incompressible neo-Hookean material 
W{Du) = /(Ai,A2) = C{\\ + A2 + A2^^A2^^ — 3) satisfies all assumptions of Theorem\^ In the case of 
a material with a linear stress-strain law, the conclusions of Theorem [7] are much more easily seen in the 
geometrically linear setting by explicitly writing down the solution v (see UP}/ )- 

The proof of Theorem [T] consists of several steps. First, we prove the existence of a sohition v for (|3.7p 
(Lemma 13.41) . Next, we show some elementary properties of v fLemma 13. 5p . which will allow us to show 
the uniqueness of v fLemma 13. 6p . Afterward we prove the remaining properties of v. This consists of the 
following steps: 

• we show that any relaxed interval has to start at i?i„ fLemma 13. 7^ : 

• we show v{Rin]T) > and v{Rin]T) < w{v' {Rin]T))Rin for some loads T fLemma l3.8p : 

• we prove that (i?i„. Rout) splits into a relaxed and a non-relaxed interval, both of them non-empty; 

• we show (|3.10p . 

Lemma 3.4. Under the assumptions of Theorem]^ there exists a minimizer v of (|3.7p . 
Proof. We start by rewriting p.7p in the form 

w, '^l'^ o ^/ r ■ friv'{r),v{r)/r) - {(p{r)v{r)y dr, 

where ip{r) = RinTin + r^^^-r^^ RoutTout is the linear interpolation between RinTi„ and RoufTout- 

We observe that {(pv)' = r (fys'^ -I- 7^'). Since f — fr for large strains by (|3.2I) and / has p-th power growth 
by (j2.8p . fr has also p-th power growth for large strains. Then the previous integral is bounded from below 

by 

C + {\v'{rr~^v'{r)+^^J -^'W^^rdr. 

The assumption p.9p is equivalent to (p'{r) > 0. Hence it is clear that the energy is bounded from below 
and that 

any minimizing sequence is bounded in W'^'P D L°°{Rin, Rout)- 

Therefore we can use direct methods of Calculus of Variations to obtain a minimizer v for this convex 
problem. □ 

In the following, we keep Tout fixed and treat Tin as a parameter, and write v(r; T) for the minimizer 
of (|3.7p with Tin = T. We also define an interval 

• — {'-^out^ ToutRout I Rin) • 

Now we prove a bound on v' , which will be useful afterward in showing the uniqueness and some properties 
of a minimizer v. 

Lemma 3.5. Under the assumptions of Theorem {1\ there exist constants 1 < Vmin < Vmax, such that for 
any T G I the minimizer v{-;T) satisfies 

Vmin < v'{r;T) <V„iax for r e {Rin, Rout)- 
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Proof. Let T G T be fixed and let v{r) :— v(r; T). We write cjr and ag for tlie corresponding radial and hoop 
stress, respectively: 

Or = difr{v' ,v/r), ag = d2fr{v' ,v/r). (3.11) 
1. Since v is minimizer of p.7|) . it satisfies Euler-Lagrange equation: 

(rcr^(r))' = ere > 0, i.e. cr'^ = - {ag - cjr) , (3.12) 

r 

where the inequality ag > follows from the definition of fr- Therefore rar{r) is non-decreasing, and 
thus we obtain 

d,U{v'{r),v{r)/T) = ar{r) > ^ar{Rrn) = — T„ > ^T,„ > ^T,^t > 0. (3.13) 

r T ^out ^out 



We see from Remark 13.31 that 

v'{r) > for r e {R,n,Rout)- 

We define 

ff(A) := 9i/(A,w(A)) for A > 1. 

This quantity represents force required to uniaxially stretch an elastic body to A times of its original 
length. It is natural to expect monotonicity of H . Under our assumptions this is true. Indeed, ()2.12p 
implies 

H\X) = (9i/(A, wiX))y = 9ii/(A, wiX)) + 5i2/(A, u;(A))«;'(A) > 

and so 

H is a, strictly increasing function. 
Using (|2.13p and the monotonicity of r(Jr{r), H{v'{r)) satisfies 

H{v\r)) < difr{v'{r),v{r)/r) = a,(r) < ToutRout / Rrn- 
So by monotonicity of H 

v'{r) < H^'^ {Tout Rout /Rin) V„iax- 

2. We want to show that cr^ > o'g in {Rim Rout)- First, let us prove that ag ^ ar everywhere. Otherwise 
let ro be such that ar{ro) — ag{ro). By differentiating (|3.1ip we obtain 

1 / , V 

V 



Cr'r = dllfr ■ v" + dufr ' 

r \ r/ 

Then using Euler-Lagrange equation p.l2|) we see that w is a solution to the second order ODE 

- {d2fr{v\v/r) - difr{v',v/r)) = di^fr{v\v/r) ■ v" + di2fr{v'.v/r) ■ - (v' - -) 
r r \ r / 

with values v'{ro) = w(ro)/ro =: k for some k. At the same time we see that V{r) :— Kr is a solution 
to the same ODE with V'{rQ) — v'{ro) and V(ro) = v{ro). Since dnfr > the ODE satisfies the 
uniqueness principle and so v — V in {Rin, Rout), a contradiction with the values of ar at Rin and 

Rout • 

Since ar ag, we have either ag > ar or ag < ar in the whole interval. In the first case p.l2|) would 
imply that ar is a, non-decreasing function of r, a contradiction with the boundary conditions for ar- 
Therefore 

ar > ag in {Rin, Rout)- (3-14) 
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3. By virtue of (P1T|) and (PTTI)) we see 

v'{r) > v{r)/r for r e {R,n,Rout)- 
Then it follows from ([^1^ and ((XT^ that 

d,fr{v'{r),v'{r)) > d,friv'{r),v{r)/r) > ^T^ut > 

and immediately 

V'{r) > Vrrnn > 1, 

where Kim depends only on /,. and RinTout / Rout ■ 

□ 

We have seen in (I3.2p that the relaxed density fr can be expressed in terms of /. As a consequence, fr 
partially inherits the strict convexity of /. We use the convexity to show uniqueness of v: 

Lemma 3.6. Under the hypotheses of Theorem\^ a minimizer v of (j3.7p is unique. 

Proof. Let u be a minimizer of (j3.7l) . By Lemma 13.51 we know that 

1 < Vrrnn < v' (r) < Vmax for r G [Rin.Rout)- 

From (|2.8|) and (|3.5p we have that (9ii/r(Ai, A2) > for Ai > 1, which together with convexity of fr in both 
variables implies the uniqueness of v' . Moreover, we know by Remark 13.11 that there exists ro € {Rim Rout) 
with nontrivial hoop stress: 

cre(ro) > 0, i.e. w(ro)/ro > w(v'(ro)). 

Since /r(Ai, A2) = /(Ai, A2) for Ai > 1, A2 > ■u'(Ai), strict convexity of / (in particular the fact ^22/ > 0) 
implies that v{ro) is uniquely determined. This together with the uniqueness of v' completes the argument. 

□ 

Lemma 3.7. Let us assume that there exists a non-empty relaxed region, i.e. there is a maximal interval 
(A,B) C {Rim Rout) such that 

v{r)/r < w{v' {r)) for r e {A, B). 
Then under the assumptions of Theorem [7] we have A — i?i„ . 

Proof. To prove that A = Rim let us assume that A > Rin- Then <tq{A) = and ag is strictly positive in a 
left neighborhood U of A. This in particular means that /(«', v/r) = fr{v' ,v/r) in U. Hence we can do all 
our computations with / instead of fr- 
We differentiate (|3.11l) in U to obtain 

'r'g = 912/ • v" + 1922/ {v' -^)= ^ (det f (y' - ^) + ^12/ • Kr) , 

where we have used ct^ = du f {v' , v / r)v" + di2f{v' ,v/r){v/ry to express v" . Now consider the limit r A. 
Since (Tg{r) — )■ 0, from the Euler-Lagrange equation rcr(, — ag — ar we know that (T^(r)r — > —ar{A) = 
-dif{v'{A),v{A)/A). Therefore we get 

<Jo{A) = (^det D'f (^v'{A) - ^) - 9i2/ • 9i/) > 0, 

where the last inequality follows from (j2.14p with Ai = v'{A) (note that by the definition of A, w{v'{A)) = 
v{A)/A). Therefore there exists no such point A and the relaxed region has to start at i?i„. □ 
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The next lemma shows that there exists a range of loads T such that v{r; T) > for r E {Rim Rout) and 
w(-;T) has a non-empty relaxed region. 

Lemma 3.8. Under the assumptions of Theorem]^ there exists a range of loads T in I such that 

< wiv'iR,mT)) andv{R,mT) > 0. (3.15) 
Proof. By Lemma [531 and (|2.10p we know that w{v' {Rin\T)) > w{Vmax), and so it is enough to show that 

< v{Rin]T) < w{V„iax)Rin 

for some range of loads T. 

Let T € X be arbitrary. As a preliminary step, we observe that T i— > v{Ri„;T) is continuous function 
for T e Z U Tout- This follows in the standard way from the uniqueness of minimizers. We also see 
that T I—)- v' [Rin]T) is a continuous function of T G I U Tout- This is an immediate consequence of the 
optimality condition difr(y' {Rin] T),v{Rin] T)/Rin) = ar{Rin) = T and monotonicity of 9i/r(Ai, A2) in the 
first variable. Then using (j2.10p we see that T 1— ?> w{v' {Rin]T)) is also a continuous function. By the same 
argument the same is true for v' {Rout',T) and w(v' (Rout'iT)) as well. 

Now we turn to the main point: the value of w(i?i„; •) at the endpoints of T. li T — Tout, the solution 
has the form w(r; Tout) — i^r for some k > 1 such that difr^K, k) — Tin — Tout, and so clearly v{Rin] Tout) = 
nRin > in this case. In fact, there is no relaxed region (and therefore no compression in the hoop direction) 
when T = Tout- 

At the other endpoint T = ToutRout/Rin we no longer have uniqueness of a minimizer for (j3.7p (see 
Remark l3.2p . Nevertheless, such minimizer is unique up to a translation, i.e. v'{-;T) is uniquely determined, 
and it has to satisfy v{Rout',T)/ Rout < w(v' {Rout]T)) . Indeed, since (rtXr)' = oe > and Rin(Jr{Rin) = 
Routes r {Rout), we have that og = 0. This immediately implies v {Rout, T) j Rout < w{v' {Rout',T)). 

Now we consider a sequence of loads Tk £ X, Tk T = ToutRout/Rin- As before we have a sequence of 
unique minimizers = v{-;Tk), where each of them has a non-empty non-relaxed region (see Remark 13. ip . 
We have shown in Lemma l3 . 71 that any relaxed region has to start at and therefore having a non- relaxed 
region implies Vk{Rout) / Rout > w{v'^.{Rout))- We also know that v{-;T) is determined only up to an additive 
constant, and so theoretically it is possible that different subsequences of {vk} are converging to different 
minimizers v{-;T). We show that this is not the case, i.e. that Vk converges to one particular minimizer 
v{-;T). Let us take a subsequence of {vk} (labeled the same) which converges to some v := v{-;T). Since 
V k {Rout) / Rout > w{v'^{Rout)), continuity of v{Rout; ■) and w{v'{Rout; ■)) implies 

v{Rout)/Rout > w{v'{Rout))- 

At the same time, we know that v doesn't have a non-relaxed region, and so v{Rout) I Rout < w{v' {Rout))- 
This shows that v{Rout) = w{v' {Rout))Rout, i-e. the limiting v{-]T) is uniquely determined. Therefore 
the whole sequence {vk} converges to this particular minimizer. This in particular shows that v{Rin', ■) is 
continuous from the left at T = ToutRout/ Rin- 
To summarize, we have shown that 

v{Rin; T) > for T = the left endpoint of I; 

v{Rin; T) < w{Vmax)Rin for T = the right endpoint of I. 

Continuity of v{Rin', •) and the fact that < w{Vmax) implies that there are some T g I such that 

< v{Rin;T) < w{Vmax)Rin- 

This completes the proof of the lemma. □ 

Corollary. Let T E T be such that (|3.15p is true. Then under the assumptions of Theorem Q] there exists 
L G {Rin, Rout) such that the interval {Rin, Rout) splits into a relaxed region {Rin,L) and a non-relaxed 
region {L, Rout)- 
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Proof. Since T satisfies p.lSp . the relaxed region is non-empty. Then by Lemma [H771 it has to be of the form 
{Rin,L) for L S [Rim Rout]- Using Remark 13. II we see that the non- relaxed region is also non-empty, and so 
necessarily L < Rout- This completes the proof. □ 

It remains to show (|3.10p . Let us compute the derivative of /i for r < L: 

h'{r) = (u))'(v'(r)) • v"{r) --( v'{r) - ^] . (3.16) 

r \ r ) 

Since r < L, the hoop stress ag = Q and together with Euler-Lagrange equation (0 — {rur)' in this case) we 
obtain 

= ae{r) = c^a/KW, u;K(r))), - = a,(r) = d^f{v'{r), w{v'{r))), 

r 

(a is a positive constant) and by differentiating 

= ai2/K(r), w{v'{r)))v"{r) + d22f{v'{r),w{v'{r))) - {w)' {v' {r)) - v"{r), 
- dnfiv'ir),wiv'{r)))v"ir) + duf{v'{r),wiv'ir))) - {w)' {v' (r)) - «"(r). 

We solve this linear system for v"{r) and (w)' {v' {r))v" {r) and substitute the result into (I3.16p . At r = L we 
have v' — v/r = v' ~ w{v'), thus h'{L) < is equivalent to 

det D^f{v'{r), w{v'{r))) - {v\r) - w{v'{r))) > 5i2/K(r), w{v'(r))) - dif{v'{r), w{v'{r))), (3.17) 

exactly matching condition (I2.14p . 

This completes the proof of Theorem [TJ 



3.2 Uniqueness of the minimizer 

In this section we show the uniqueness of the minimizer for the relaxed problem (|3.1I) . 

First, it is easy to check that uq defined in (|3.6p is a minimizer for the relaxed energy Eq (since the 
functional is convex, any solution of the Euler-Lagrange equations is a minimizer). Since Eq is not strictly 
convex, the uniqueness of the minimizer is not clear. It is however true: we now show that uq is (up to 
an additive constant) the only minimizer of the relaxed problem. Indeed, suppose that there exists another 
minimizer ui of the relaxed problem. By convex duality for Eq we have: 

min Eo{u) ~ max / —D{(T)dx', 

diver— 0, (T.n— TatcJfi 

where I?(cr) := supp-gjjaxs (cr, F) — Wr{F) is the convex conjugate of Wr- Since uq is a minimizer of Eq, we 
know that the maximum on the RHS is attained for ctq = ^^(^"o)- From the definition of the convex 
conjugate D we see that {aQ,Dui) — Wr{Dui) < D{a-Q), and after integration Eq{ui) > J^~D{a) dx' . 
Since the deformation ui is a minimizer of Eq as well, we obtain an equality in the last relation. Hence 
{ao,Dui) — Wr{Dui) = D{aQ) a.e. in fl and consequently ctq = ^gp-{Dui) a.e. in fl. It follows that 

DuQ = Dui (3.18) 

at points where Wr is strictly convex at Duq (i.e. at points where the eigenvalues Ai > A2 of {Duq Duq)^/^ 
satisfy Ai > 1 and A2 > u;(Ai)). 

We have proved in Theorem [T] that the hoop stress is tensile exactly in the non-relaxed region 

nN := {x:L< \x\ < Rout} (3.19) 
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with Rin < L < Rout- By fJ/f = £7 \ ri^r we denote its complement, i.e. the relaxed region. Since the both 
stresses are tensile in fijv, as a consequence of p.2p we have Wr{DuQ) = W{Duq). 
Using the strict convexity of Wr in the non-relaxed region n jv we see that 

uo{x) = ui{x) + C (3.20) 

ioi X € flN- We may assume without loss of generality that C — 0. In the relaxed region we need to 
replace p.lSp by 

Duq{x) ■ n{x) = Dui{x) ■ n{x), x G ^Ir 

where n{x) = j|j is a unit vector in the radial direction. This concludes the proof of uniqueness, because 
integrating the last relation in the radial direction extends validity of ()3.20|) to the whole set ft. 

4 The 2D result 

In this section we show the matching upper and lower bounds in the simplified Kirchhoff-Love setting. After 
stating the main result of this section we prove the upper bound by superimposing wrinkles on the solution 
uq of the relaxed problem. We use Lemma HTT] to create simple wrinkles whereas Lemma UH] provides a tool 
to create a family of wrinkles with changing wavelength near the free boundary. In the rest of the section we 
prove a matching lower bound. Using Lemmas 14.41 14.61 and 14.71 we show that if Eh{uh) is close to fo, then 
also Duh has to be close to Duq. As a consequence we obtain a bound on the out-of-plane displacement Uh,3,- 
By interpolation we show the smallness of Duh.s, which allows us to project Uh to the plane without changing 
its energy too much. The final ingredient is a comparison of uq with the projection of Uh (Lemma 14. 8p . 

For Theorem[2]we will assume that the lower bound in (|2.8p holds with p = 2 rather than just 1 < p < 2. 
The stronger assumption p ~ 2 is only required for the second half of the proof of the lower bound (e.g. for 
the interpolation), whereas the first half of the proof of the lower bound (especially Lemmas 14.41 14.61 14.71 
and the Poincare inequality) requires only 1 < p < 2. Since the real purpose of the two-dimensional result 
is to lay the ground for the proof of the three-dimensional case, many of the preparatory lemmas are proved 
in the more general setting 1 < p < 2. 

Let us now state the main result of this section: 

Theorem 2. Let us assume (|2.8p with p — 2 and all hypotheses of Theorem]^ Then there exist constants 
< Ci < C2 independent of h such that 

£0 + Cih < min Eh{u) <£o + C2h, (4.1) 

u 

where £0 is the minimal value of the relaxed problem l|3.ip . 

Remark. Theorem\^ doesn't necessarily require all hypotheses of Theorem]^ In fact, some of them can be 
replaced by assumptions on the solution uq of the relaxed problem (which are consequences of Theorem\^. 

4.1 The upper bound in the two-dimensional setting 

To obtain the upper bound, we must construct a test function Uh for any (small) h > Q with energy 

Eh{uh) < £0 + Ch, 

where the constant C is independent of h. A naive approach would be to superimpose a "single family of 
wrinkles" (with a well-chosen period independent of r, and a well-chosen amplitude that depends on r), 
c.f. [TU] (see also [TT] for a similar calculation). The energy associated with this Uh has the expected scaling 
away from r = L (the edge of the wrinkled region). However the membrane and bending energies are both 
singular at r = L; as a result, the total energy (after integration) is too large, of order £q -f 0{h\ \ogh\). The 
singularity in the bending term can be avoided by introducing a boundary layer, however we have not found 
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a similar way of avoiding the singularity in the membrane term (which is associated with stretching in the 
radial direction). 

To get a linear correction a more complicated construction seems necessary, using a "cascade of wrinkles" 
rather than a "single family of wrinkles." In other words, the period of the wrinkling changes repeatedly as 
one approaches the edge of the wrinkled region. Constructions of this type have been used in other settings, 
for example in studies of compressed thin film blisters [H [12] . 

Recall that the solution of the relaxed problem has compressive hoop strain when r < L. The essential 
purpose of the wrinkling is to avoid this compressive hoop strain by out-of-plane buckling. In the following 
lemmas, we write e*(r) for the compressive hoop strain to be avoided. Up to a factor of 27rr, this amounts 
to "the amount of arclength to be wasted by wrinkling" along the image of the circle of radius r. 

Lemma 4.1 (see Lemma 2 in [4 ). For every e, > there exists a smooth C°° planar curve 7(e*) = 
(71 (e*)(0) 72 (£*)(*)) : M — > with properties 

\da\ = 1 + e„ dai > 0, j{-t) = --fit) 



7(i + 27r) =7(i) + 







and satisfying the bound 



l7i - ^1 + \dai - 1| + \dtai\ < Ce„ I72I + 19*721 + 19**721 < CeV\ 



|4.7i| < C, \d,,-f2\ < Ce-'/\ \dl-f2\ < Ce, 



-3/2 



where C does not depend on e,. Moreover, the hound is sharp (in terms of the scaling in e<,) for small values 
of e^. 



Idea of the proof (see Lemma 2 in for more detail): We construct 7 by reparametrizing the curves 

7(t) = 



t 

As'mt 



so that |9t7| = 1 + £*, where A is chosen such that J^^ \dt"/{t)\ dt = 27r(l + By considering the small- 

1/2 

limit we obtain A — p(e*)e* (p being a smooth function on [0,1]). This leads easily to the desired 
estimates. □ 



We would like to use Lemma 14.11 to superimpose wrinkles on top of the planar deformation uq obtained 
from the solution v to the relaxed problem (13. 7|) . Though this naive construction does not achieve optimal 
energy scaling, it can be modified (using Lemma 14. 2p to obtain a construction with optimal energy scaling. 
Therefore it makes sense to analyze the naive construction, to understand why it fails and also to motivate 
the successful construction. 

The naive construction is based on Lemma l4.1l and proceeds as follows. After we obtain the parameter 
e*(r) (amount of wastage of arclength) from v, we determine the right period of wrinkles to obtain the 
optimal scaling. A bit of calculation reveals that the stretching and bending terms are of the same order 
when the number of wrinkles is of order 

k := [h-'/% 

where the brackets denote an integer part. From v we obtain the amount of wastage of arclength as 

v[r)/r 

for r G {Rim L). Following the proposed idea we define a solution Uh in (using radial coordinates r, 9) by 

- r / ^» 71 (£*(?■)) (^6*) , .l2{e,{r)){ke) 

Uh{r,e):=:v[r)r+ ' ^ '^"^ ^ 9 + v{r) ' ^ ^ 63 (4.2) 
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and as uq elsewhere. From p.lOp wc know that ^ L — r (up to a factor) for < L — r « 1. Since 
72 ~ £* "^^^ by Lemma |4.1[ the contribution of drrUh,^ to the bending energy is divergent: 



I \drrUh,3\^rdrr^h^k-^ f ({L - r)-^'"^^^ r dr = [ {L-r)-^rdr. 

It is clear that the bending energy is of order h in the region Rin < r < L ~ h. Therefore a boundary layer 
in the region L — h < r < L would solve this issue provided the integral over this layer of the new \drrUh,3\'^ 
is at most of order h^^. 

Now let us try to compute the contribution from the stretching energy near the transition from the 
relaxed to the non-relaxed state. In the region Rin < r < L — S (for a fixed 6 > 0) the term drUh.3 from the 
membrane energy satisfies 



L-S nL-S 

\drUh,3\^rdr^ / k-^\v'{r)j2f + Hr)d,, rdr. 

By Lemma 14.11 the first term in the parentheses is of order e* whereas the second term is of order e~^. 
Therefore, near the free boundary (where e, << 1) the second term is dominant and we obtain 

\drUh.3\^rdr^ \ k-h:^rdr^hi\og{5)-log{L-R,n)). 

We see that by setting S — h {oi any power of h) we would obtain energy scaling h\ logft,|. 

We now begin discussion of the successful construction, which uses a "cascade of wrinkles" near r = L. 
The main tool is Lemma 14.21 Whereas Lemma |4 . 1 1 involved wrinkled curves. Lemma 14.21 involves wrinkled 
strips, in which the length scale of wrinkling doubles from one side to the other. It provides the basic building 
block for our cascade of wrinkles. 

Lemma 4.2. Let B = (0, 1) x (0, w) with < /, w < 1, F : (0, /) — > K, and e be a positive function on (0, 1) 
satisfying \e'\ < c, \e"\ < c, and l/c < e < cl for some c > 0. Then there exists a smooth deformation 5'(s,t) 
defined on B and w-periodic in the t variable such that for any t S (0, w) the following holds: 

^i{s,t) = F{s), se{0,l) 

i<f2, *3) (s, t) = w-f{e{s))i-), s e (0, //4) 

w 

(vI/2, *3) is, t) = |7 (e(s)) (^^) , . e (3/4 1, 1) 
4'(s,0) = (F(s),0,0), ■^{s,w)^{F{s),w,Q), s e (0,1), 



and 

|9.*2| V l^.^-aP < Cw^ri, \dt^\^l + e{s), \D^'f\^<C{l-^w^ + w-~^l), 
where C depends just on c and 7 is the curve defined in Lemma \4-l\ 

Proof. To prove this lemma we just need to define 5*2, ^'3 such that the required estimates are true. The 
idea of the construction is very similar to the proof of Lemma 14.11 To simplify the notation we first assume 
w = 1. 

First, let us fix < £ < 1 and a e [0, 1]. We consider a planar curve: 

J : t ^ {t,A[{l - a)sin(27rt) + asin (47rt)]) , 
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where A = A{e, a) is such that the length of 7([0, w]) is exactly (1 + s). More specifically, we define 
A{t) = -^ I |7'(r)|dr = [ J 1 + (27rA)2 [(1 - a) cos(27rr) + 2a cos(47rT)]^ dr 

with A such that A(l) = 1 (since A is strictly increasing function of A, there exists a unique such A). 
Considering the small-e limit, we obtain 

A = /)(e,a)£^/^, 
where p G C°°{[0, 1] x [0, 1]). Consequently we have that 

dA 



Using 7 we define a new reparametrized curve 

T{s,a,t) =70 A-i(f). 

This curve obviously satisfies |r'(t)| — 1 + e and 

ri(e, a,t) =t - spi{e,a,t), T2{e,a,t) = e^^'^p2{s, a,t), 

where pi,p2 € C°°([0, 1] x [0, 1] x M). From there we get estimates 

\d,T,\ + \dtT,\<C, \d„Ti\<Ce, 
l^eFil + IdeaTil + Id^tTil < C, |5„„ri| + |a,*ri| + IdttTil < Ce, 
la^^a^SJTal < Ce^/^-fe k,l,n> 0. 

Now we are ready to define the map ^. We set 

:=F(s), *2(s,i) :=ri(e(s),<^(s),t), *3(s,i) := r2(e(s), .^(s), i), 

where ^ is a. smooth increasing function on (0, /) satisfying 

./)(s)=0 sG(0,Z/4), 
0<(/>(s)<l se (1/4,3/4 0, 
<^(s) = l s€ {3/4 1,1), 



and (j)' <3/l. Then 



5«^'(s,t) = (F'(s),aerie' + a«ri<^',aer2e' + a„r2.^')^ 
at^'(s,t) = (o,atri,5tr2), 



and using previous estimates together with e « / we obtain desired bounds 

\dt^\ = 1 + e(s), |5«*2|' + |5«vl/3|2 < Cl-\ 

To finish the proof in the case w = 1, we get the estimates on D"^^ in the same way as for the first derivatives 
of*: 

iD^'fil < C 

|5„*2|<c;, |a,t*2|<c, \dss'^2\ < cr^ 
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It remains to show the lemma for general w, i.e. we need to define on B — (0, 1) x (0, w). Let ^ come 
from the proof of this lemma for the case w = 1 is defined on (0, 1) x (0, 1)). Then we simply set 

^{s,t) = (^F{s),w'^2{s,t/w),w'^3{s,t/w)j . 

Now it is an easy calculation to show that ^ satisfies all the estimates. □ 

Remark. Later we will use Lemma \4.S\ to create a test function for the upper bound for Ef^. The definition 
of the function F from Lemma \4.^ will be based on v (the solution of the relaxed problem p.7p ]. and 1 + e 
will be the corresponding natural width. Let M{s) = F'{s)ei ® ei + {1 + e(s))e2 ® e2 be a 3 x 2 matrix. 
Since 1 + e(s) is the natural width, we have that DW{M(s)) = a{s)ei ei for some scalar function a(s). 
By definition we know M(s)ii = ds'i'i{s). Then the boundedness of D^W (see ()2.8p ) implies 

[ W{D^)- [ W{M{s))< C [ 1(9,*, at*)-(F'(.s)ei®ei+at^' (8)62)1^ 

JB JB JB 

= / \ds^2\' + \ds^3\^ 
JB 

where we used that = 1 + e (see Lemma \4-.2^ and the rotational invariance of W . Using Lemma \4-S\ 

we obtain 

[ WiD'f) + h^\DH\'^ - [ WiMis)) <C {w^ + [w^r^ + w-H^]) . 

JB JB 

As already mentioned, we will later set F :— v, the solution of the relaxed problem, and e := e*, the excess 
arclength, so that J W{M) is the energy of the relaxed solution. This remark will be then used to compare 
the elastic energy of the constructed deformation with the energy of the relaxed solution (which is Sq). 

Remark. Since * is periodic in t, we can assume it is defined in an infinite strip (0, 1) x M. 

Remark 4.3. In Lemma \4-.2\ we have estimated the size of , not only D^^^ (in fact, the third component 
of 4* was the most oscillatory, and so it is larger than other two). We will use this fact later in the proof of 
the upper bound in the general three-dimensional setting. 

By p. 101) and smoothness of v we can choose a small 6 > and constants < C2 < ci such that 
— ci < 9r£* ^ —C2 in the interval (L — S, L). We define a deformation Uh by changing Uh (defined in (j4.2l) ) 
in the region L — S < r < L. The idea is to create a cascade of wrinkles by superimposing wrinkles coming 
from Lemma 14.21 in smaller and smaller rectangles as we approach r = L. Wc define (for a non-negative 
integer n): 

/„ := (L - 54"", L - (54-("+i)), 
a„:=L-54-", |/„|, u;„ := |^2-", 

and for r G In set: 

Uh{r, e) = v{r)f + + u(r)*3£3 
where ^ comes from Lemma [4.21 applied to the rectangle /„ x (0,w„) with 

F(s) w(s - a„), e(s) := £*(s - a„), w := u;„, I := L- 

We do this construction in the region Un := [j^^gln x (0,27r) with N — — l/21og2(/i) {N is chosen such 
that 2^^ = h~^). Since the doubling of a period defined in Lemma [4.21 happens strictly inside the given 
interval, the first and second derivatives of Uh are continuous at each a„ (i.e. there are no jumps in the first 
derivative of between /„ and 

To finish we just need to define Uh in a region close to r — L. Consider the strip (L — 64^^ ,L) (observe 
that it includes In)- By Lemma [4.11 the amplitude of Uh,3 at gn is of order WN£*'^^^iaN) ~ b^l'^h?!'^ . The 
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length of the interval In is 5 (4~^ — 4-(^+i)^ = 3/4 . j/j, "We multiply Uh by a smooth cut-off function in 
In to bring the out-of-plane displacement of Uh to zero at a^+i- Since the length of the interval In is of 
order h and the value of Duh is bounded in that interval, the membrane energy is bounded by Ch in this 
interval. The second derivative of the new deformation in this region is at most of order h~^/^, and so the 
bending energy in this region is less than Ch'^ * {h-^^^f Ch^. OveraU we obtain: 

N 

Eh{uh) - Eoiuo) <Ch + Cj2'2"k [wl + {l;;^wl + uj-Hl)] 

n=0 

<Ch{l + h^2^'^) < Ch. 

4.2 The lower bound in the two-dimensional setting 

In this section we want to prove the lower bound 

min Eh(u)>£o+ch (4.3) 

for some c > independent of h. Our argument uses the convexity of the relaxed problem; we shall have to 
work a lot because the relaxed problem is not strictly convex. We will proceed by contradiction, assuming 
there is a deformation u with energy very close to the energy £q of the relaxed solution. After obtaining a 
bound on the out-of-plane displacement we use interpolation to show smallness of Du^. This allows us to 
project the deformation u into the x-y plane without altering its energy too much (i.e. we obtain a planar 
deformation with energy close to £0). Finally, we conclude the proof by showing that it is not possible to 
have a planar deformations with energy close to Eq. 

The results in the first part of this section will be useful also later for the proof of the three-dimensional 
case; therefore we prove them assuming only 1 < p < 2 in (|2.8p . On the other hand, it is convenient to 
assume p — 2 in the interpolation argument used later in this section, and so from that point on we assume 
p = 2. 

We define 

for some 1 < p < 2. We observe that the function gp is monotone, convex, and C^. Since it also satisfies 
gp{2t) < 4(7p(i), convexity of gp implies 

gp{a + b)<2{gp{a)+gp{b)). (4.5) 

The proof of the lower bound is divided into six steps: 
Step 1: To proceed by contradiction, we assume that for any small < 6 < 1 there exists a sequence of 
functions Uh s.t. 

Ehiuh) < £0 + Sh. 

This is equivalent to 

Eh{uh) - Eoiuo) < Sh, (4.6) 
hence using Wr < W and definition of Eo we immediately obtain 

WriDuh) - WriDuo) dx' + B{uh - uo) = Eo{uh) - Eo{uo) < Sh. 

Since uo is the minimizer of the relaxed energy Eo, it has to satisfy Euler-Lagrange equation 

/ DWriDuo) : Dip Ax' + B(^) = 



n 
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for any test function £ VF-'^'P(r2, R'^). Since DWr{DuQ) is bounded, we can easily take test functions in 
W^-'^'-'^(f2, K^). Using the relation for <fi := uo — Uh yields 

[ Wr{Duh) - Wr{Duo) - DWriDuo) : {Duh - Duq) dx' < Sh. (4.7) 
Jn 

Step 2: We would like to obtain a pointwise lower bound on the integrand of the last relation. 

Lemma 4.4. Assume (|2.8p holds. Let F e M"^^^ be an orthogonal matrix with singular values Ai > 1, 
•^2 < Ai, and n be the right singular vector corresponding to Ai (i.e. F'^ Fn — Xln). If there exist k > and 
an open neighborhood U o/(Ai,A2) such that 

D'^fr{oi,(T2) > KCi (g) ei for {ai,a2) £ U, (4.8) 

then there exists Cg > (depending only on Ai, A2, k,U, and growth of W) such that for any G G R'^^^ 

Wr{G) - Wr{F) - DWr{F) : (G ~ F) > co5p (|(G - F) ■ n\) . (4.9) 

Proof. We first prove the statement for large strains (i.e. if G is large). Let ai > a2 are the singular values 
of G. We observe that the coercivity of W{G) (see (12. Sp ) implies that Wr{G) has also p-th power growth 
for large matrices. Indeed, if G has only tensile stress, we have W{G) = Wr{G), i.e. Wr{G) has the same 
growth as W. Otherwise, we know from the ordered force inequality (j3.4p that 

/r(ai,a2)>/.(ai-l,l)>c(|Gr-l), 

where the last inequality follows from the coercivity of W and from the fact that strains (Ai — 1,1) produce 
only tensile stresses (and so = Wr in this case). Thus the LHS of (|4.9p grows at least like \G\p. Since the 
RHS has at most such growth, the conclusion follows. 
It remains to prove the statement in the case 

\G\ < M 

for some M. First, we observe that in this case \{G — F) ■ and \{G — F) ■ n\P are comparable, and so we 
can replace gp in (j4.9p with a quadratic function, i.e. we need to show 

Wr{G) - Wr{F) - DWriF) : {G ~ F) > ca\{G - F) ■ (4.10) 

for any \G\ < M (possibly with a different cq than in (|4.9p ). We start by computing DWr{F). Since Wr is 
rotationally invariant, we can assume without loss of generality that 

Then a simple calculation reveals that 

fai 0\ 

DWriF) = a2 , 

\o Oj 

where ai = 9i/r(Ai,A2) and a2 — d2fr{M, ^2)- We observe that ai > 0, a2 > 0, and that ^2 = iff 
A2 < w{Xi). We rewrite (|i?TUl) : 

WriG) - Wr{F) + DWriF) : F > DWriF) : G + coliG ~ F) ■ (4.11) 

- a,Gn + 02622 + Co ((Gn - Ai)' + Gj, + G^^) 

= (ai - 2coAi) Gn + a2G22 + cq {Gj^ + Gl^ + G^i) + cqA? 
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We choose cq > small enough so that ai — 2coAi > 0. The LHS of the inequality depends only on 
singular values of G. Hence, we can prove (j4.1ip by maximizing the RHS among all matrices G with given 
singular values ai > ui- We give the argument assuming that ol\ — 2coAi > a2 (the proof in the case 
ai — 2coAi < 0-2 is analogous). 

To prove (|4.1ip we will use the following lemma due to von Neumann (see, e.g, [H]): 

Lemma 4.5. If A, B are nxn matrices with singular values 

(Tl > • • • > Cr„, pi > ■ ■ ■ > Pn 

respectively, then 

n 

\tT{AB)\ < ^arPr- 

r=l 

We shall apply the lemma to find the maximal possible value of the RHS (I4.1ip among all matrices G 
with singular values ci > (T2- First we set A := (G|0) (the 3x3 matrix with first two columns identical with 
G and the third column equal 0) and B = diag(ai — 2coAi, ck2, 0). Then the lemma gives 

(ai — 2coAi) Gil + a2G22 < (ai — 2coAi) cti + a2CT2- 

For A := G^G and B = diag(l, 0) the lemma implies 

G\i + G21 + G31 < tJi. 

Together we see that the RHS of (|4.11l) is at most [ai — 2coAi)cri + €^2(^2 + cqo-J + cqA^. To see that this 

/ai 0\ 

bound is optimal, we use Go := (72 as G. We got that the RHS of (|4.1ip is maximal for the choice 

Vo 0/ 

G = Go, and so we need to prove (|4.10|) only for G — Gq for any ci > (T2 > 0, i.e. to show 

/r('Ti,(72)-/.(Ai,A2)-i?/.(Ai,A2)(ai-Ai,a2-A2) >co((Ti -Aif. (4.12) 
Using Taylor's expansion of /,■ at the point (Ai, A2), and assumptions (|2.8p and (|4.8p we see that 

/r(^7i,a2)-/.(Ai,A2)-i?/,.(Ai,A2)((Ti-Ai,a2-A2) > " '\^'/r(0(f^-e),<^-Ode 

"'(Ai,A2) 

> cik((Ti - Ai)^, 

where the last inequality follows from the fact that D^fr > nei (g) ei in a non-trivial part of the segment 
between (Ai,A2) and (o'i,cr2) (remember that G is bounded), and D^fr > otherwise. This completes the 
proof of the lemma since we showed that ()4.12|) holds with cq = cik. □ 

Lemma 4.6. Assume (|2.8p holds. Let F G M'^^^ be an orthogonal matrix with singular values 1 < A2 < 
Ai < K. If there exist k > and an open neighborhood U of (Ai, A2) such that 

D^fr{cri,a2) > kI for (di, 0-2) e 

then there exists cq > (depending only on Xi, X2, k,U, K , and growth of W ) such that for any G £ R^^^ 

Wr{G) - Wr{F) - DWriF) : {G - F) > c^g^ (|G - F|) . (4.13) 

Proof. The idea of the proof is simple and resembles proof of the previous lemma. If G is large, we get the 
statement the same way as in the previous lemma. 
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Otherwise we can assume \G\ < M for some (possibly large) M. For any such G the part of the segment 
connecting F and G (i.e. F + t{G — F) for t £ (0, 1)) which belongs to U will be at least e part of the whole 
segment. The LHS of (|4.13p can be written as 

/ {D^Wr{F + t{G - F)){G -F),G^ F) (1 - t) dt, 
Jo 

and the integral is at least k(1 — t)\G — F\'^ along the non-trivial part of the segment and non- negative 
everywhere else. Therefore (I4.13P follows. □ 

Before proceeding to step 3, we need one more lemma, that is similar in character to the preceding ones 
but involves W instead of Wr- 

Lemma 4.7. Assume (|2.8|) holds. Let F G R^^^ be an orthogonal matrix with singular values Ai > 1, 
A2 < w{Xi), and n be the right singular vector corresponding to Ai. // there exists k > and an open 
neighborhood U o/(Ai,A2) such that 

£'^/r(o'i,o'2) > KCi (g) ei for (cti,CT2) e U, 

then there exists cq > such that for any G G K^^^ 

W{G) - W{Fo) - DW{Fo) : (G - F) > coQp (dist(G, ^0(3)Fo)) , (4.14) 

where Fqu = Fn and Wr{F) = W{Fo). 

Proof. Let H E M.^^"^ be such that Hn = Gn and W{H) is minimal among all such H. We observe 

W{H) > WriH), W{Fo) = WriFo). 

Arguing as we did in the proof of Lemma |4]4] (i.e. using rotational invariance of W and Wr to compute DW 
and DWr) we also have 

DWiFo) -.(G-F)^ DWr{Fo) : {H - Fq). 

We claim that 

WiG) - W{H) > cgpi dist(G, SO{3)H)) 

follows from (|2.8|) . We give a proof of this fact in the Appendix (see Lemma [A. 31 and Remark I A. 4^ . 
Now Lemma 14.41 and the previous inequality imply that the LHS of (|4.14p is at least 

[W{G) - WiH)] + [Wr{H) - Wr{Fo) - DWriFo) : {H - Fq)] > 

CO {gp (dist(G, SO{3)H) + gp {\{H - Fo) ■ n\)) (4.15) 

(cq is a generic constant, i.e. it can change from line to line). 

We claim that in the second term on the RHS of (|4.15p we have 

\iH ~ Fo) ■ n\ > cdist(i7, 50(3)^1)). 

To prove it, we first observe that without loss of generality we can assume n — (1,0). Then using (jA.7P we 
see that 

H = {vi\v2) and Fo = (wi|w2), 

where fi _L D2, |f2| = ff(|wi|) and wi _L W2, \w2\ — g{\wi\). There clearly exists a rotation R G 5*0(3) such 
that and V2\\Rw2, and the vectors have the same orientation; thus 

dist(i7, 50(3)Fo)2 < \vi - i?U.iP + \V2 - i?W2p = (It'll - klD' + (|f2| - \W2\)^ 

= (Kl - \w,\r + igi\v,\) - gi\w^\)f < c(K| - \w^\r < c\iH - Fo) • nf. 
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Hence the RHS of (I4.15P satisfies 



CO {gp (dist(G, SO{2,)H) + 9p (I - ^o) • n\)) 

> Co [gp (dist(G, SOiS)H) + g^ (dist(i/, SO{3)Fo))] > c[,5p(dist(G, 50(3)Fo)), 



where we used incquahty ()4.5p . This completes the proof of Lemma 14.71 □ 



We continue the proof of the lower bound (|4.3|) . 
Step 3: It is clear that there exists p > I such that the region where the smaller singular value of Duq is 
at least p is non-empty. We will denote this region Np and its complement in il as Rp. Lemmas 14.41 and 14.61 
applied to (|4.7p imply 

gp {\Duh{x') - Duq{x')\) dx' + [ gp {\{Duh{x') - Duq{x') ■ n{x')\) dx' < CSh, 



where C > depends only on uq, Wr, and choice of p. We now use the Poincare inequality adjusted to our 
setting (Lemma IA.6p . We obtain 

9p{\uhix') - Ch - uq{x')\) Ax' < [ gpi\Du{x') - Duoix')\) dx' < CSh 
for some Ch- 

We would like to extend the previous estimate into the whole f7. We fix a direction 9 and define (using 
radial coordinates) 

f{t) Uh{t, e) - uo(i, 9), t e {Rtn,Rout)- 
Let us call M the radius of the boundary between Np and Rp, and let 

K:= r"^\p{\rmdt+ r°^\p{\f{t)-cu\)dt. 

jRi„ JM 

Then by the Poincare inequality applied to / on {Rm, Rout) we get 

gp{\f{t)-cg\)dt<C r^^' gp{\rmdt < CK 



R,r^ J Ri 



for some ce- We also have 



{Rout - M)gp{\ch ~ ce\) = I gp{\ch -~ ce|) dt 



M 



and so 



<C\ gp{\f{t) - Ch\) dt + / gp{\f(t) - ce\) dt\<CK 



Ro^t ( fRa.,t fRa^t \ 

9p{\f{t) - Ch\) dt<c[ / gp{\f{t) - ce\) dt + / gp{\ce - c,,\)dt < CK. (4.16) 



R,: 



Finally, by integrating (|4.16p in 6 we obtain 

gpQuhix') -CH- uo{x')\) dx' < CSh. (4.17) 
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step 4: The next step in the proof is the interpolation between ||u/i^3||L2(n) and \\D^Uh,3\\L2{n)- For that 
reason we assume p — 2 (instead of a more general 1 < p < 2) for the rest of this section. (Since the previous 
lemmas will be used later (see Section[S]), we proved them assuming only 1 < p < 2.) 
When p = 2, (|iT7)) reads 

\\uh{x') -Ch- uo{x')\\l2(^n) dx' < CSh. 

Without loss of generality we can assume Ch = 0, since our problem is translation invariant; in particular we 
have 

IMlhin) < CSh. (4.18) 

Since 

Eh{uh) - h^\\D'^Uh^'i\\\2i^Q^) > mmEa{u) = Sq, 

we have that 

\\DW,,\\hin)<S/h. (4.19) 
Interpolating between (|4.18p and (|4.19p we obtain 

WDu^Mnn) < CS. (4.20) 

Step 5: We want to estimate 

' \WiDuh)-WiDul^)\ dx', 



where ujj^ = {uh,i,Uh,2, 0) is the projection of u/,. into x-y plane and A{x') := Duh{x') — Duj^{x'). By (|2.7p 
the previous integral is equal to 

' \f{a,{x'),a2{x')) ~ f{Xi{x'),X2{x'))\ dx' , 

where ai{x'),a2{x') and Ai(x'), A2(a;') are the singular values of Duh{x') and Duj^{x'), respectively. Since 
the singular values of a matrix are Lipschitz functions of the corresponding matrix (see, e.g., Corollary 8.6.2 
in [12]), we have 

\f{a,{x'),a2{x')) - f{X,{x'),X2{x'))\ = \Df{^{x')){a^{x') - X,{x'),a2{x') - X^ix'))] < C|I)/(C(x'))||A(x')|, 

(4.21) 

where ^{x') is a point on the segment connecting {ai{x'),a2{x')) and (Ai(x'), A2(x')). By (12.81) Df{l, 1) = 
and D^ f < C, and so Df{C) < C(|C| + 1)- Using quadratic growth of / ( (^75)) with p 2) we obtain 

\Df{a^')t < Cm^')\' + 1) < C'ifi^ix')) + 1) < C'{f{<j,{x'),a2{x')) + f{X^{x'),X2{x')) + 1), (4.22) 

where the last inequality trivially follows from the convexity of /. Integrating (|4.2ip and using (|4.22p together 
with Holder inequality we get 

\f{a^{x'),a2{x'))~f{X^{x'),X2{x'))\ dx' < C ( [ \Df{ax'))\'dx'X^'(f \A{x')\' dx'^ 



< C 

Therefore by (|4.20p and using S < l,h < 1 we have that 



.1/2 

J^^W{Duh) + W{Dul^)dx' + l\ \\DuhM\LHn)- 



( [ \W{Duh) -W{Dul^)\ dx'^ ^^(2 [ W{Duh)dx'+ [ \W{Duh) - W{Dul^)\ dx' + l) S 
\Jn J \ Jn Jn J 

<C6(^£o+Sh + l + J \W{Duh) - W{Dui^)\ dx'^ < C'S (^1 + ^ \ W{Duh) - W{Dul^)\ dx' 
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It follows easily that 

/ \W{Duh) - W{Dul^)\ dx' < CS^/^, 
Jn 

and so 

\Eh{uh) - Eh{ul^)\ < Cb"!^ + h^Wuh.zWl^^a) < C(^'^' + Sh) < C5^/\ (4.23) 
Equations (|4.6p and (|4.23p together imply 

E{ul')-Eo{uo)<CS'/^ 

where E(v) — W{Dv) + B(v) (v is just an in-plane deformation, so there is no bending term present). 

Step 6: Observe that the last relation does not depend on thickness h anymore. In fact, to finish the proof 
we just need to show that the minimum of the energy Eh over in-plane deformations has to be strictly larger 
than the minimum of the relaxed energy Eq'- 

Lemma 4.8. Under the assumptions of Theorem\^ we have 

min Eh{u) > min Eq(u) = £q. 

Proof. Let us assume the contrary, i.e. for any small S > there exists a function u : — > such that 

/ W{Du) dx' + B{u) <£o + 5= [ Wr{Duo) dx' + B{uo) + S. 
Jn Jn 

The plan is to obtain a contradiction by showing that the areas of u(Ofl) and uo{i^R) should be very 
similar using one argument and at the same time very different for another reason (here ^Ir is the relaxed 
region introduced near p. 191) ). First, using Euler-Lagrange equation for uq we can replace the boundary 
term B(u — uq) by the gradient term: 

[ WiDu) - WriDuo) - DWriDuo) ■ {Du - Duq) dx' < S. 
Jn 

Since W > Wr and Wr is convex, the integrand in the last relation is non-negative a.e. Therefore the last 
relation remains true if we integrate over the relaxed region VIr instead of the whole domain fl. To proceed, 
we would like to find a matrix Fq such that Wr{Duo) = W{Fo) (i.e. we want to relax compressive stresses 
in Duq if they are present). We know that Duq has compressive stresses in the hoop direction (and tensile 
in the radial direction), which means that Fq and Duq coincide in the radial direction and are different in 
the hoop direction, i.e. 

Fo{x')n{x') = {Duo{x'))n{x') and F„{x')n-^{x') = cn-^{x') 

with some c > {Duon^)n^ . Moreover, it can be easily seen that in this case DWr{Duo) = DW{Fo). We 
rewrite the previous inequality to obtain 

/ W{Du{x')) - W{Fo{x')) - DW{Fo{x')) ■ {Du{x') - Duo{x')) dx' < CS. 
Jnu 

Using Lemma 14.71 with p = 2 we see that 

/ dist^(i:'w, 5*0(3)^0) da;' < C5. 
Ja.R 

Whence 

\u{nR)\~ f detFoda;' 

Jnrc 



<[ \detDu~detFo\ dx' < C (S + S^/A <CS^/^. (4.24) 
Jnn ^ ' 
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Further we see that 

detFo > detZ^Mo. (4.25) 

Indeed, Fq and Duq are "same" in the radial direction (i.e. Fqu = {Duo)n), whereas due to compression in 

the hoop direction in the relaxed solution uq we have {Fon^)n^ > ((Duo)n^)n-^ . Hence we get (|4.25p by 
taking product of the two previous relations. Finally integrating (|4.25|) we obtain 

/ detFoda;'> / det Dug dx' = \uoiflR)\. (4.26) 

To finish the proof we want to show 

\u{nR)\ < \uo{nR)\ + CS^^\ (4.27) 
Then by combining (|4.24p and (|4.27p we get 



/ detFoda;' < |uo(f^i?)| + CJ^/^^ 



contradicting (I4.26P since S > can be arbitrary small. 
To show (|4.27p we set e := (5^/^ and define 

M := u{flR) n{xeR^ : dist(x, uo(^7i^)) > e} . 

Then 

\u{nR)\ <\{x£R'^ : dist{x,uo{nR,)) <e}\ + \M\< |uo(f^fl)| + Ce + |X| . (4.28) 
It is enough to estimate the size of A4. Arguing as before (c.f. (I4.17p ). we know that 

II" ~ '"o|li2(Oj,,R2) < C6, 

whence 

e^\u~\M)\< I \u-uo\^dx' <\\u~uo\\l2,nj,R2)<C6. (4.29) 

By virtue of (lOO)) 

\M\^ detDudx'= / detFodx'+ / det - det Fq dx' 

Ju-^(M) Ju-^(M) Ju-^(M) 

<C\u-^{M)\+ I IdetDw-detFol da;' < C(5e"2 + C(5^/2 < C(5^/^ 

Using this estimate in (|4.28p we obtain (|4.26p . □ 
This completes the proof of the matching lower and upper bound in the two-dimensional setting (|4.ip . 



5 The 3D result 

In this section we will prove the scaling law for the minimum of the elastic energy in the nonlinear three- 
dimensional setting. As in the previous section, we need to show an upper and a lower bound. As usual in 
problems of this type, the upper bound is an easy consequence of the upper bound for the 2D setting. The 
main goal is therefore to show the lower bound in this more general setting. 
As explained in Section [2] we consider a nonlinear 31? energy 

Er{u):^\( WMDu)dx 
Jnx{o,h) 
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instead of the reduced 2D energy (|2.2p . The boundary term in the 3D setting is defined as 

(«) ^ / .^dS-^ f u{x) ■ dS. 

The main resuh of this section is 

Theorem 3. Under the hypothesis of Theorem]^ there exist constants < Ci < C2 independent of h such 
that 

£o + Cih< min E^" (u) + B^^ (m) < £0 + C2 /i. 

uew^-p{nx{-h/2,h/2)) 

5.1 The upper bound in the three-dimensional setting 

The construction of a test function uf^^ in the three-dimensional setting is based on the test function 
defined in the previous section to show the upper bound in the two-dimensional setting. Following proposed 
Kirchhoff-Love ansatz, the normal to the mid-surface remains straight and normal to the mid-surface after 
deformation. Therefore, we just need to find how much should each of these normals stretch. It follows from 
the definition of W (see (12.51) ') that for any x' E there exists a factor a{x') such that the vector a{x')i'{x') 
satisfies 

W{Duh{x')) = W3D{Duhix')\aixyix')), 

where iy{x') is the unit normal to Uh(ri) at Uh{x'). We observe that a{x') is bounded and \Da\ < C\D'^uii\- 
We define the solution u'f^^ as 

u^{x) := Uh{x') + xy, ■ a{x')v{x') 

and compute 

Duf{x) = {Duh\ av) -I- • {D{av)\{)). 

Then 



- EhWh) <T [ Wsoi {Duh\ai^)+X3 ■ (D(az/)|0) ) - WMDuh\ai^) dx 

lnx{-h/2,h/2) 



< \ ( • DW:iD{Duh\av) : {D{av)\Q) + x\ ■ C\D[ptv)\^ dx 

IQ.x(-h/2,h/2) 



<Ch^ I \D{av)\'^ dx <Ch'^ [ a{x'f\Di^f + \Dafdx' 
Jn Jn 



where C depends on HD^WsdHloo. We know that \D^uii\^ is bounded by Ch (see Remark H3)) . and 

thus using boundedness of Duh and 

/ \Dv\^dx' <C [ \D^Uh\^dx' 
Jn Jn 

we obtain h^ a^\Dv\'^ dx' < Ch. A similar estimate is true for the second term: 

h^ I \Da\'^ dx' < Ch^ [ \D^Uh\^ dx' < Ch. 
Jn Jn 

Together we have obtained 

En^r) < Eh{uh) +Ch<£o + Ch. 
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5.2 The lower bound in the three-dimensional setting 

Our goal is to show a lower bound similar to (|4.3p : 

un .o™}'\ /o , ^h"" i^) + {u) > f o + ch. (5.1) 

To first approximation, the proof of (jS.ip consists of slicing our domain, applying the two-dimensional lower 
bound on each slice, then patching them together. But two new features require changes in the argument. 
First, the energy density W in 2D was derived from the 31? energy density W^d assuming the missing third 
component is optimal. Therefore we need to estimate how much the third component of Du differs from the 
optimal one, and how the optimal vector depends on the first two components. As a second feature, where 
2D used interpolation we need to proceed differently by using a rigidity theorem. 

From now on let /i > be fixed and consider a function u e x (— /i/2, h/2)). Then for any fixed 

CC3 € {—h/2, h/2) we define 

u{x') :— u{x' , X3) 

for x' (z fl. Since uq is the minimizer of the relaxed energy Eq we have that 

£0 = Eoiuo) < Eo(u) < / W{Du) dx' + B{u) 
Jn 

and we set 



R := [ W{Du) dx' + B{u) -£q>Q. 
Jn 



Our initial goal is to show that if R is small then Du is close to some function (which is derived from 
Duq). We use Euler-Lagrange equation to replace boundary term B by DWr and obtain 

/ W{Du) - Wr{DuQ) ~ DWr{DuQ){Du - Duq) dx' = R. 
Jn 

Using notation of Lemma 14.71 this can be rewritten as 

/ WiDu) - W{Fo) - DW{Fa){Du - Duq) dx' = R, 
Jn 

where Fo{x) ■ n{x) — Duo{x) ■ n{x) and W{Fq{x)) — Wr{Duo{x)). We apply Lemma W7J\ to get 

, ( diBt{Du{x'), SO{i)Fo{x'))) dx' < CR. (5.2) 



n 

Now we go back to the 3D body. We see that 

rh/2 



El'^iu) + Bl'^iu) I I WMDu{x)) - W{Du{x)P) dx- 

" J -h/2 Jn 




W{Du{x)P) dx' + B{u{-,xs,)) - £0 

Du(x') I 



dx3 =: h + h 



'\ 0^ 

where P = I 1 1 . Let ^[x) £ g^^j^ ^hat 



,0 Oj 



W3DiDu{x)P\C{x)) = W{Du{x)P). 



28 



Then Lemma IA.2I implies 



h 

and by virtue of (|5.2 



h>^ I ' I 9p{\Du{x) - {Du{x)P\ax))\)dx, (5.3) 



-h/2 Jn 



h>Y / gp{dist{Du{x)P,SO{3)Fo{x')))dx. 

" J-h/2 Jfl 

det {Fo{x')\V) > and ' 



c /■''/^ 

'h/2 Jfl 

We want to extend ^0(2;') G R'^^^ into a 3 x 3 matrix. To do that, we find a vector V E M.'^ such that 

W3DiFoix')\V) = mmW^oiFoix'm. 
Observe that the choice of V is unique. We define 

M{x') {Fo{x')\V). (5.4) 

Lemma lA.51 then impUes 

dist((i:'u(a;)P|C(a;)),50(3)M(a;')) < C dist{Dil{x)P, SO{3)Fo{x')), 

and so 

/2>^/ f 9p{dmt{{DiL{x)P\C{x)),SO{3)M{x')),)dx. (5.5) 



h/2Jn 



By adding (j5.3p and (|5.5p . and by using (|4.5p (the triangle inequality for 5^ with factor 2) we obtain 



El^{u) + Bl^{u) - £0 > X / 9p( dist(L>t2(a;), S'0(3)A/(x'))) dx. 



h 



nx{-h/2M/2) 



Moreover, from the previous analysis (see (j4.17p ) we know that there exists a function t of 3:3 alone such 
that 

i / / gp{\u{x)-T{xs)-uo{x')\)dx'dxs<C{El^{u) + Bl''{u)-£o). 
" J~h/2 Jn 

By adding those two inequalities we obtain 

Jin,u)<C{El°{u)+Bl^{il)~£ol 

where for U d we define 
J(C/, u) :== 

gp{dist{Du{x),SO{3)M{x'))) dx + l- f gp {\u{x) ~ t{x3) ~ uo{x')\) dx. 

''■ Jux(-h/2M/2) " J(7x(-ft/2,/i/2) 

We will obtain the lower bound from the following important lemma: 
Lemma 5.1. Let h > be sufficiently small and let vq be such that (ro,ro + 2h) C {Rin,L). Then 

where ,4(a,/3) is an annulus with radii a < (3 and rj = 7?(ro) > is a decreasing function of r^. 
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The proof of the lemma will be given in Section 15.31 The desired lower bound is an easy consequence. 
Indeed, for any i?i„ < vq < L — 2h we know by Lemma IS.ll that 



J{A{ro,rQ + h),u) >?7(ro)/i^ 
Adding such inequalities for ro = Rin + kh such that Ri„ < ro < {Rin + L)/2 — 2h we obtain that 



K K 



J{A{R,n,L),u) > J2 v{roik))h^ > V i{R^n + L)/2) J2h^>Ch 

for some C > 0, where K = {Rin — L)/2h — 2 and we used monotonicity of rj. 

Besides proving Lemma (see the next section), the proof of the lower bound (|5.ip is finished. 

5.3 Proof of Lemma 15.11 

Proof. Let us first sketch the idea of the proof. We assume u has small energy (i.e. J{u) is small), and want 
to compare it with uq. We take a collection of neighboring cubes in the hoop direction, each cube 

with side h. Using a rigidity theorem we show that ii is almost constant on each cube and doesn't change 
much between the cubes. 

Since u has small energy, we prove that often Du is larger than Duq in the hoop direction. After 
integration in the hoop direction we obtain that u — uq can not be small in most of the cubes. On the other 
hand, u having small energy implies that u — uq has to be small in the sense, contradicting the previous 
fact. 

For better understanding we split our proof into several steps. 
Step 1: Let us consider a part of A{ro, tq + h) with the length approximately 2h^/^. We set 

k := [/1-1/2] 

and for a given 9 g (0, 2tt) we define 2k cubes in the radial coordinates 

= {tq, ro + h)x{e + ia, 9 + {i + l)cr) x {-h/2, h/2), 
where i = 1, ... ,2k and a = h/rQ. We denote by Q the union of those cubes and set 

J ■.^hJiQ,u) + \Q\h^ 



= / gp{dist{Du{x),SO{3)M{x'))) dx+ / gp i\u{x) - T^xa) - uo{x')\) dx + \Q\h^ . 

Jq Jq. 

Step 2: We claim that J > Ch?-^ for some positive C. To prove the claim, we shall suppose that 

J = e/l^■^ (5.6) 

and give a lower bound for e. Since M is defined in terms of uq (see (|5.4p ). we have \DM{x')\ < C and 
consequently 

gp {dist{Du{x), S0{3)M,)) dx < C [ gp {dist{Du{x), S0{3)M{x'))) dx + C\Q^\h^, 

where Mi := M{xi) with Xi being a point in Qi (e.g. center of Qi). Using the rigidity estimates of [T] (see 
also [Tl]) we obtain a rotation Ri on each cube Qi such that 



gp i\Du{x) - RiM,\) dx < CJ (5.7) 
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and 



2fe-l 



J2 gp{\R^-R^+l\) < CJh 

i=l 

Using convexity of gp and Jensen's inequality we get 

\Ra - Ra+l3+l\ \ ^ Y.t=a 9p {\Ri - Ri+l\) 



< < CH 



/3 



'0.5 



(5.8) 



Since the RHS in the previous relation is smaU, so is the argument of gp on the LHS (in particular, it is 
smaller than 1), and so we can take square root of both sides to obtain 



We also know 



i?, e S0{3), \M,\ < C, \n\ = 1, \M, - Mj\ < Ckh, \n{x) - n{y)\ < Ckh. 

Therefore we can choose one rotation i?* (e.g. among Ri), a matrix Af* (among Mi) and a unit vector 
such that for any x' 

\R^M,n{x')^ - R^M^n:t\ < C {6 + kh) . (5.9) 
Step 3: For j e {1, . . . , fc} we have 



u{x) — t(x3) — uo{x') dx — / u{x) — t{x3) — uo{x') dx 



Qj+k 



{Du{x) - Duo{x')) n{x'y if j{x')dx 



j [{Du{x) ~ R^M,) n{x')^ + {R^M,n{x')^ - RM*ni) + 

{RM*ni - Duoix')n{x')^)] ipj{x') dx, (5.10) 



where Dj :— [jjl^j Qi and < ipj{x') < h is a. weight coming from the integration (more precisely, ipj is 
linear in Qj going from to h, ^Pj(x') = h on Qj+i, . . . ,Qj+k-i, and then decays linearly from /i to in 
Qj-i^n)- We point out that used in (|5.10p means a unit vector in the orthoradial (hoop) direction. The 
first two parts of the last integral can be directly estimated from (|5.7I) and (|5.9I) using convexity of gp and 
Jensen's inequality as in (|5. 



{Du{x) ~ R^M,) n{x')-^(p.i{x') dx 



< 



and 



{R,M,n{x')^ - R^M^.ni) ip-j{x') dx 



<C\Dj\ {Sh + kh^ 



(5.11) 



(5.12) 



Step 4: Now we show a lower bound for the remaining term Jj^ (^R^,M^n^ — Duo{x)n{x')'^) ipj{x') dx using 

the fact that Duo is "smaller" in the hoop (n^) direction than R^M^,. Arguing as in the proof of ()4.25p we 
see that for any unit planar vector x we have 



|M*X| - \Duo{x')nix')^\ > k > 
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where k depends monotonically only on the radial position r = \x'\ (and approaches as ro — > L). Therefore 
(using Ko ■■= K{ro + h)) 



(R^M^n^ — Duo{x)n{x')^) fjix') dx 



> 



Ri,M^n^ / ipj{x')dx 

JDi 



Duo{x')n{x')^ipj{x') dx 



>\M^n^\ [ ipj{x')dx- [ \Dun{x')n{x')-^\(pj{x') 
= / [iM^n^l ~ \Duo{x')n{x')'^\) ipj{x') dx 

JDj 

> / Koipj{x') dx > Koh\Dj\/2. 



dx 



Step 5: Using ()5.1ip and ()5.12p together with the last relation we see from (|5.10p that 



u{x) ~ 'T'ix'i) — Uo{x') dx 



u{x) ~ T(x-i) — wo(x') da;J > 

\D.j\h (ko/2 -8-kh)- ChJ^'^\Dj\^'^. (5.13) 



From dlH) and 5 = VCJkh^ « ^/jJF^ < y^, we see that ko/2 - 6 - kh > /vo/4 > for e < Ckq and 
small h. 

To finish the argument we sum (|5.13p over j = 1, . . . , to obtain 



Ck\Q\h{Ko/2-5-kh)< / \u{x) - t{x3) ^ uo{x)\dx + ChkJ^^^\Q\^/^, 



(5.14) 



where we have used that \Q\ and \Dj\ are comparable. Using the convexity of gp and Jensen's inequality we 
have 

" \uix) - Tixs) - uoix)\dx < ^i/2|Q|i/2, 
and so after plugging the values oi S, \Q\, and k into ()5.14p we see that 



Hence we obtain 



J > CkIH^-^ and J(Q, u) > Cnlh^-^. 
Finally we cover the annulus yl(ro,ro + h) with approximately 27r7'o/i^^/^ distinct copies of Q to obtain 

J(^(ro,ro + h),u) > ?7(ro)/i^, 

where rj^ro) = CvoKq. □ 

6 Discussion 

We have identified the scaling law for the minimum of the energy. We have achieved this by (i) constructing 
a family of functions with low energy; (ii) proving a lower bound on the energy with the same scaling as in 
(i). 
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In the construction for the upper bound we introduced a cascade of wrinkles to obtain the hnear scahng 
in thickness. As we approach the boundary between relaxed and non-relaxed region, the amplitude of the 
out-of-plane deformation vanishes, and so does the bending energy. The decay of the amplitude balances 
the increase in the bending energy due to the increasingly fine scale wrinkles. 

The "naive construction" discussed at the beginning of Section U has no cascade, but its energy is 
So + Ch\ log h\ rather than £q + Ch. It is natural to ask whether something similar to our cascade is required 
to get the optimal scaling. 

This paper has focused on the energy scaling law. It is natural to ask how the energy is distributed more 
locally; for example, is the optimal distribution (with respect to radius) similar to that of the construction 
giving our upper bound? It is equally natural to ask what the minimizer looks like; for example, must the 
amplitude and wavelength of wrinkling at radius r resemble these of our construction? 

It seems worth noting that while we have repeatedly used the Euler-Lagrange equation for the relaxed 
problem (to characterize the relaxed solution) , we never used the Euler-Lagrange equation from the original 
problem. In related studies, such as |20) . minimizers are known to have special properties. We expect the 
same to be true in the present setting, but the analysis of minimizers will require new techniques. 



Appendix 

In this section we prove several lemmas which were used previously. Lemma [A.2l and Lemma lA.SI were used in 



the proof of the lower bound in the three-dimensional setting f Section 1 5. 2 1) . Poincare inequality (Lemma lA.GP 
was used in the proof of the lower bound in the two-dimensional setting f Section 14. 2p . Lemma lA.31 and 
Remark I A . 41 were used in the proof of Lemma l477l and Lemma fA.ll was used in the definition of W (see (|2.6p ). 

Lemma A.l. Let W3D be a stored energy function of an isotropic elastic material with 

where J := det(F), C := F'^ F , and 

h := J-'" tr{C), h ^ {itr{C)f - tr(C^)) . 



// g satisfies 



^(/l,/2,J)>0, ||(/i,/2,J)>0 (A.l) 



for Ii > 3, /2 > 3, and J > 0, then for any M G M^^^ we have 

mmW3D{M\0= min W3d(M|0, 

where {M\S,) denotes a 3 x 3 matrix with first two columns identical with M and the third column ^, and 
M _L ^ means the columns of M are perpendicular to ^. //, moreover 

|^(/i,/2,J) + |^(/i,/2,J)>0 (A.2) 
0/1 oh 

for Ii > 3, I2 > 3, and J > 0, the minimum of WsoiMl^^) is attained only if M J- S,- 

Proof. Isotropy of the material implies rotational invariance of the energy density W3D, and so without loss 
of generality we can assume that 

% 0' 

Af = I A2 I . (A.3) 
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Therefore we want to minimize 

/Ai 
W^D A2 2/ 
\0 

among all possible x, y, z. A simple calculation reveals that 

J = A1A22;, 

h = (A1A2Z)-'/' (A? + Xl + x^ + y^ + z^) , (A.4) 
h = (AiAaz)--*/^ {XlXl + Xfz' + Xjz' + Xjy' + Xjx^) . 

Using the AM-GM inequality we see that /i > 3 and /2 > 3 even if x = y = 0. Fixing z and varying x and 
y the value of J stays constant whereas /i and I2 are increasing functions of x^ and y^. Therefore (jA.l[) 
implies that W^d has its minimal value (for any z fixed) if x = y = 0. In the case (jA.2[) the conclusion 
follows by the same reasoning. □ 

Lemma A. 2. Let W^d satisfy (|2.4I) and (12.61) . Then for any M £ R'^^^ there exists a constant C > such 
that 

WMM\v) - W{M) = WMM\v) - WMM\0 > Cgp i\v - ^\) = Cgpi\{M\v) - (M|C)|), 
where ^ — argmin H/3£)(Af| •). 

Proof. If |u — ^1 is large, the conclusion follows from the growth condition (|2.4p . Let us therefore assume 
that |u — ^1 < i^. As in the proof of the previous lemma we may assume M satisfies (|A.3|) . We write 

v^ix,y,zf, e=(0,0,Zf, 

where f can be written in this form due to the previous lemma. We have 

W3d{M\v) - W3d{M\^) = {W3d{M\v) - W3d{M\w)) + iWMM\w) - WMMIO) , 

where 

w := {0,0, zf. 

We estimate 



(A.5) 



W^D {M\v)- W^D {M\w)= g{J, h,h)- g{J, h,h) 

- V.g • (0, J-2/3(a.2 + y2)^ J~^l%X\y^ + Xlx^)) > C{x^ + y^), 

where J, Ii , I2 are defined in (|A.4|) and 

h = (A1A2Z)-'/' {Xl +Xl + z^) , 
h - (A1A2Z)-'/' {XjXl + Xlz' + Xlz^) . 

We also have 

WMM\w) ~ W3u{M\i) = f3D{Xl,X2,z) - /3d(Ai, A2,Z) 

= 933/315 (Ai, A2, C)(^ - Z)V2 > C(z - Zf, 

where we have used that (93/3li(Ai, A2, Z) — and D'^fsD > (see (12.41) ). Adding (jA.5l) and (|A.6|) we obtain 
the desired inequality. □ 



(A.6) 



Lemma A. 3. Let the density function W^d satisfy (j2.4p and p.6p . For any unit vectorui £ K and another 
vector vi £KP, \vi \ > 1, we define 

F argmin {WsoiF) : F e R^''^,Fui = i;i,detF > O} . 
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Then there exists an orthonormal basis ui,U2,U3 and orthogonal vectors wi,W2,W3 such that 

Fui^Vi, 1 = 1,2,3, 
\v2\ = \v3\ =g{\vi\), 



(A.7) 



where g{t) is a Lipschitz continuous function. Moreover, there exist < Ci — Ci(ui) < C2 — C2(wi) such 
that 

Cigp (dist(G, 50(3)F)) < WMG) - W^F) < C2 dist2(G, 50(3)F), (A.8) 



for any G satisfying Gui — vi, det G > (function gp was defined in (j4.4p ). 

Remark A. 4. As a consequence of (IA.8I) we obtain a similar condition for W - for a given unit vector 
n G and another vector m e R"^ we have 

Cigp{dist{G,S0{3)F)) < W{G) - W{F) < G2dist^(G, 5*0(3)^), 

when Gn = Fn = m and W{F) — miii{W (H) : Hn = m}. 

Proof of Lemma \A.3[ We first prove (jA.7[) . Since W3D is frame- indifferent and isotropic, we can assume 
WLOG that ui — ei and vi = AiCi, Ai > 1. Tlien we want to find all matrices F that minimize 

miii{W3D{F) : Fei — Aiei}. 

Let F be such a matrix. We see from Lemma lA.ll that the first column of F and the third column of F 
are perpendicular. Since Fei = AiCi, this means that F12, = 0. Since the second and third column are 
interchangeable, the first and second column are also perpendicular, i.e. F12 — 0. By the same lemma we 
also have that the second and third column of F are perpendicular. Therefore, up to a rotation which fixes 
the first column, F is diagonal: 

/Ai ' 
F= A2 
V Aa^ 

We claim that A2 = A3. Indeed, strict convexity of fsD implies 

/3d('^1, A2, A3) + /3d(Ai, A3, A2) 



A2 + A3 A2 + A3 , 



W^3d(^^) = /3I?(Al,A2,A3) = 

> /3d(Ai, ^^^y^, ^^^^) = WMF') 
(where F'ei = AiCi), unless A2 = A3. Moreover, since 

A2 = argmin /3D(Ai,t,t), 
t>o 

it follows that A2 is a Lipschitz continuous function of Ai. 
We observe that Ai > A2. Indeed we will show that 

/3d(Ai,Ai,Ai) </3Z3(Ai,t,0 (A.9) 

for any t > Ai. We compute 

J(Ai,t,i) ^ Ait^ > A? = J(Ai,Ai,Ai), 
h{Xut,t) = Xt^h-^^' + 2A-'/'t2/3 > 3 ^ ^^(^^^ 
l2{Xi,t,t) - X-^^'t^/' + 2A?/V2/3 > 3 = /2(Ai, Ai, Ai), 



^See 1 12. 6I I for definition of J,Ii, and l2- 
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where the two latter inequalities follow from the AM-GM inequality. Therefore for i > Ai we have that 
Ii{Xi,t,t) > /i(Ai,Ai, Ai), l2{Xi,t,t) > /2(Ai, Ai, Ai), J{Xi,t,t) > J(Ai,Ai,Ai), and so (|A.9p follows 
from (ED). 

To prove (jA.Sp . we fix the matrix F (therefore Ai > A2 are also fixed). It is sufficient to show the lower 
bound and upper bound only for G for which W3d(G) — WsoiF) is small: 

Cidist'^{G,S0{3)F) < W3d{G) - WsoiF) < C2dist^ {G, SO (3) F) (A.IO) 



(the upper bound and the lower bound for "large" G follow from the growth assumptions on f^o ^ see (12.4^ ). 

We start with the lower bound. As before, we can assume F is diagonal with entries Ai, A2, A3 = A2. Let 
CTi > cr2 > 173 > be singular values of G. We first see that 

Ai = Gil = tr(Gei 18) ei) < cti • 1 = cri, 

where the inequality follows from von Neumann's Lemma (Lemma 14. 5p . and Gn denotes the upper left entry 
of matrix G. Since WsniG) depends only on the singular values of G, we will try to estimate both sides 
of (jA.Sp in terms of the singular values of F and G. 
We start with estimating LHS of (|A.10p . We know 

dist^iG, SO (3) F) ^ mm ti {{G - RFf {G ~ RF)) 

i?G SO (3) 

= tr (G^G) + tr (F^F) - 2 max tr (RFG^) . 

-R(z5'0(3) 

We claim that maxfi^so{3) tr [RFG^) is equal to the sum of singular values of FG^ . To show this, we use 
the singular value decomposition FG^ = UDV^ , where U,V G 5*0(3) and D is a diagonal matrix (with the 
singular values of FG'^ as diagonal entries). Then tr {RFG'^) = tr {RUDV^) = ii {V'^RUD) = tr {QD). 
Since Q is a rotation, we have that tr {QD) < tr (D), and the maximum is attained (for Q — I). 
The sum of singular values of FG^ is by definition equal to 

tr((TO^GF^)'/') tr( (G^G)'/V 

We write F — (Ai — A2)ei ® ei + A2/ and obtain 



tr 



((G^G) f) = (Ai - A2) (s/Wg) + A2tr (\/g^) = (Ai - A2) {Vg^^ + A2 {cji +02 + 0^). 



We know that vG^Y? has eigenvalues cti, 172, ca, and so tr(G"^G) = (7\ + (72 ^ '^3- Therefore we see that 
dist2(G, 50(3)^) = {af + cri + al) + (A? + A^ + A^) - 2A2(cri + da + ^3) - 2(Ai - A2)a 



(A.ll) 

= (ai - Ai)^ + (^2 - A2)" + (as - A2)^ + 2(Ai - A2)(ai - a). 



where a denotes 



a = (VG^)ii = (VG^^ei, ei). (A.12) 

We claim that ^ 

a>^. (A.13) 
0-1 

Indeed, we write G — U DV, where U,V G 50(3) and Z? is a diagonal matrix with entries ai,i = 1,2,3. We 
define a unit vector 

X = {xi,X2,X3) := Vei. 

By virtue of (IA.12P we get 

3 

a = (\/G^ei, ei) = {V'^DVei^ei) = {DVei, Vei) = ^ xlui. 
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We also know that Gei ~ XiCi, and so 



Aici = UDVei = U aiXjej^ 



In particular, since U e S0{3), we see that the norm of the vector on the RHS is Ai, i.e. 

3 



Y.a^af = >^. (A.14) 



i=l 



To summarize, we have 



a = ^^x^(Ti, (A. 15) 

4=1 

^x2 = land ^x^af ^A?. (A.16) 

1=1 i=l 

To find the lower bound for a, we simply optimize J2^=i ^l^i assuming (jA.16p . Using method of Lagrange 
multipliers we observe that for X]i=i ^I'^i to be minimal one Xi has to be zero. We see that X]i=i ^I'^i is 
equal to: 

Ai + (72(73 -r n -^1 + ^1^3 -f „ ^1 + ^1^2 .f „ 

11 Xi — 0, it X2 = 0, if Xs — 0. 

(72 +(73 (7l+(73 (7l+(72 

Using convention (7i > (72 > (73 and (jA.14[) we see that (7i > Ai > (73. Since 

Af+(7l(73 Af+(72(73 2^ 2 , ^1 +'^l'^3 , Af +(7l(72 2 ^ 2 

S ^ (73 and < A]^ — "'1 5 

(7l + (73 (72 + (73 (7l + (73 (7l + CT2 

the minimum of ((XT5|) is equal to ^^7^^- Finally, we observe that a\ > implies 

Aj +(71(73 ^ A| 
(7l + (73 ~ (7l 

We have proved (|A.13p . 

Now we are ready to finish the proof of the lower bound. By virtue of (jA.lip : 

dist2(G, 50(3)^^) < ((71 - Ai)2 + (^2 - \2f + ((73 - A2)2 + 2(Ai - A2) ( ai - ^ 

For the middle term in (|A.10|) we have: 

WsoiG) - W-iD{F) = /3d((7i,(72,(73) - /3Li(Al, A2, A2) 

3 

> ^9i/3_D(Ai, A2, A2)((7.j - Xi) + C I ((71,(72, CT3) - (Ai, A2, A2)|^ 

i=l 

3 

= ai./3Z5(Ai, A2, A2)((7i - Ai) + C^((7, - A,)', 

2=1 

where we have used strict convexity of /3D (see \2A\ ) and 92/3d(Ai, A2, A2) — d^fsuiXi, X2, A2) = (a con- 
sequence of the definition of A2). Therefore, it remains to show 

2(Ai - A2) (ai < CdihD{Xi,X2, A2) ((71 - Ai). (A.17) 
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We observe that since cti > Ai, we have ^cti — < 2{ai — Ai). Smce Ai,A2 are fixed and dif^o > 0, 

(|A.17p immediately follows (with constant C depending on Ai). 

We now turn to the proof of the upper bound. This is easy, since we have already done all the work. We 
know that 

dist2(G, SOiS)F) - (ai - Ai)2 + (as - Aa)^ + (^3 - A2)' + 2(Ai - A2)((Ti - a) 
> (ai - Xif + ia2 - A2)2 + (^3 - X2f + 2(Ai - A2)((7i - Ai), 

where we have used a < Ai. This is true by (|A.15I) and (IA.16P : 

3 / 3 \ / 3 \ 

^a;,(x,a,)< 5]a;n ^a^^^H = Ai. 



a — 



We also know 

W3d{G) - W-iD^F) = /3D(CTl,0'2,cr3) - /3D (Ai , A2, A2) 
3 

- X! ^^-^31) (Al, A2,A2) (cTj - \i) + C I ((71,(72,^3) - (Ai, A2, A2)|^ 

2=1 

3 

= ai/3Z5(Ai, A2, A2)(ai - Ai) + C^((T, - A,f , 

2=1 

where we have used boundedness of D^fsjj (see (|2.4p ). To finish the proof, it remains to observe that 

ai/3D(Ai, A2, A2)(<Ti - Ai) < C(Ai - A2)((Ti - Ai) 
holds trivially (with Ai > A2 being fixed). □ 

Lemma A.5. Let F,G e R^"^^. Let ^ € satisfies ^ -L F, det(i^|^) > 0, and \^\ = /(F), and similarly let 
e satisfies C _L G, det(G|(') > 0, \(^\ = 1{G), where 1{A) is a Lipschitz continuous function of singular 
values of A. Then there exists constant C such that 

dist((F|C),50(3)(G|C)) < Gdist(F,50(3)G). 

Proof. Without loss of generality we assume columns of F lie in the x-y plane, i.e. F31 = F32 = 0. Let 
R e 50(3) be such that 

dist(F, S'0(3)G) = \F- RG\. (A.18) 

We show that R can be chosen such that columns of RG lie in the x-y plane as well. We know that (see the 
proof of Lemma IA.3P : 

dist2(F, 50(3)G) = tT{F^F) + tr(G^G) - 2tr ((FG^GF'^)^/^) . 

Since Fji = F32 — 0, the RHS in the last relation does not change if we replace F by its first two rows. The 
last term is then equal to tr ((G^G)i/2(F'^F)i/2) , and we obtain 

dist2(F,50(3)G) = ||(G^G)i/2 _ (F^F)i/2||2. 

Now let S e 5*0(3) be such that (5G)3i = (6*0)32 = 0. We want to compute 

dist(F, 50(2)5*0), 
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where we treat F and SG as 2 x 2 matrices (since both their third rows vanish). Fohowing the previous 
reasoning we obtain 

dist^ {F,S0{2)SG) = WiiSGfSGy/^ - {F^Fy/^\\^ = \\{G^Gy/^ - (F^Fy/^W^ 

We have shown that dist^(F, 50(2)5'G) = dist^(F, 50(3)0), i.e. that R can be chosen such that RG lies 
in the xy plane. 
We have 

dist^((f^|e),50(3)(G|C)) < \\{F\0-R{G\0\\' 

^\\F- RG\\'^ + - i?CP = dist2(F, 50(3)0) + |C - i?CP- 

Since F and RG lie in the x-y plane, both ^ and R^ are perpendicular to this plane. It is straightforward but 
tedious to show that in fact ^ and i?C have the same orientation. Then we just use the fact that |^| = 1{F) 
and Id = 1{G) together with Lipschitz continuity of I to obtain 

\£,-RC\ = \l{F)-l{RG)\ < C\F-RG\ = Odist(F, 50(3)0). 

We are done since (|A.18[) and the previous relation imply 

dist2((i^|0, 50(3)(0|C)) < |(F|0 - i?(0|C)|' - |F - RGf + \^ - RCf < Cdist''{F, 50(3)). 

□ 

Lemma A. 6 (Poincare inequality). Let gp be as in (|4.4I) with 1 < p < 2. Then there exists a constant 
C{U,p) such that for every v E W^'^{U) there exists a constant v and 



p{\v-v\)dx < C / 5p(|Vw|)dx. 
Ju 

Proof. We first observe that since gp{t) < ^min{tP,t'^) and gp is convex, there exists O such that 

9p{s + t) <C{sP + t^), for every > 0. (A.19) 
In the proof we will use the following truncation result proved in [14) : 

Proposition (Truncation). Suppose U C M" is a bounded Lipschitz domain. Then there exists a constant 
G{U,p) with the following property: For each v G W^'P{U) and each A > 0, there exists V : U ^M. such that 

(i) llVl^lUo. < OA 

(ii) \{xeU: v{x) ^ V{x)}\ < § 4,^^|v.(.)|>A} I Vt'r dx, 

(iii) l|V^;-V^/|^^,(^) <0/^^^^^l^^(^)l>,j|Vz;|^'dx. 



Let us denote 



K -= I gp{\Vv\)dx. 

u 

l.oo 



By the proposition with A = 1 there exists V G W such that |VV^| < O and 



iVw- V^H^p < O / \Vv\Pdx<CK. 

{\Vv\>l} 



The standard Poincare inequality implies 

[ \V -V\^dx <C [ \VV\^dx <C i [ \VV\^dx+ [ \VV\^ dx 



lU Ju \J {v^V} J {v=V} 

<C\{v^V}\+C 



I \VV\'^dx<GK. (A.20) 

J{v=V} 
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We also get 

/ \V-v-a\Pdx<C [ \VV -Vv\Pdx <CK. (A.21) 
Ju Ju 

Using (|XT9| . (|X20l) . and (|X2T|) we obtain 

J gp {\v - {a + V)\) dx<c(^J \V - dx + J \V - v - a\P dx^ < CK. 

□ 
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